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■ Abstract 

We study the Dirac sea in the presence of external chiral and scalar/pseudoscalar 
, potentials. In preparation, a method is developed for calculating the advanced and 

retarded Green's functions in an expansion around the light cone. For this, we first 
expand all Feynman diagrams and then explicitly sum up the perturbation series. The 
light-cone expansion expresses the Green's functions as an infinite sum of line integrals 
over the external potential and its partial derivatives. 

The Dirac sea is decomposed into a causal and a non-causal contribution. The 
causal contribution has a light-cone expansion which is closely related to the light-cone 
ON ' expansion of the Green's functions; it describes the singular behavior of the Dirac sea 

in terms of nested line integrals along the light cone. The non-causal contribution, 
on the other hand, is, to every order in perturbation theory, a smooth function in 
position space. 
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Qh. 1 Introduction 

<D 

In relativistic quantum mechanics with interaction, the fermionic wave functions are 
solutions of a Dirac equation of the form 

X" 

^ , (i$ + B-m) * = , (1.1) 

where B is composed of the classical bosonic potentials. According to the common con- 
ception, the Dirac sea of the system is built up of all the negative-energy solutions of the 
Dirac equation. We can describe it with the so-called fermionic projector P Q. On the 
non-rigorous level of this introduction, the fermionic projector is given by the formal sum 
of the projectors on all these solutions; i.e. 

P(x,y) for = Uy ^ a (x)%ti) , (1.2) 

where the index 'o' runs over all the quantum numbers of the negative-energy states. 
We want to analyze how the fermionic projector depends on the bosonic potentials in 
fll.ip . According to the decomposition Ql.2j ) into the individual states, this dependence 
can be regarded as a collective effect of all the fermions of the Dirac sea moving in the 
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external potential B. Following Dirac's original concept that the completely filled Dirac 
sea describes the vacuum, we can also say that we are interested in how the fermionic 
vacuum is influenced by the bosonic fields. Our aim is to describe this physical effect in a 
detailed and explicit way. 

It turns out that the dependence of the fermionic projector on the external potential has 
a complicated non-local structure. In order to simplify the problem, we shall study P(x, y) 
in an expansion about the light cone, which is called light-cone expansion. The light cone 
around a space-time point x consists of all points which can be reached from x with a light 
ray. In flat Minkowski space, which we will consider here, the light cone is given by all pairs 
of points [x, y) of Minkowski space whose Lorentzian distance (y — x) 2 = (y — x)j (y — x)- 7 
is zero. In the simplest case of a smooth function A(x,y), the light-cone expansion is just 
an expansion in powers of (y — x) 2 , i.e. a representation of the form 

oo 

A(x, y ) = Y,(y~ x ) 2j M x >y) (1-3) 

3=0 

with smooth functions Aj(x,y). Since the expansion parameter (y — x) 2 vanishes on the 
light cone, the coefficients Aj{x,y) give approximations of A(x,y) in a neighborhood of 
the light cone (i.e., Ao(x,y) coincides with A(x,y) on the light cone, A\(x,y) gives the 
first order behavior of A(x,y) for pairs (x,y) which are close to the light cone, etc.). The 
important point is that the Aj(x,y) are approximations of A(x,y) even for points x,y 
which are far apart. We only need that the pair (x, y) is close to the light cone, which 
is an unbounded hypersurface in IR 4 x IR 4 . In this sense, the light-cone expansion is a 
non-local expansion. The major advantage over local approximation techniques (like e.g. 
Taylor expansions in the space-time coordinates) is that the light-cone expansion gives a 
much more detailed description of the fermionic projector in position space. Furthermore, 
since the light cone is the boundary of the domain of causal dependence, all the effects 
related to the causality of the Dirac equation occur near the light cone. Thus the light- 
cone expansion describes the fermionic projector precisely in the region which is most 
interesting physically. In this paper, we will develop an efficient method for performing 
the light-cone expansion of the fermionic projector. 

After this simplified and very qualitative introduction, we briefly discuss the difficulties 
and methods of the more detailed study. First of all, it is not obvious how to characterize 
the "negative-energy solutions" of the Dirac equation in the case with general interaction. 
In other words, one problem is to find the right quantum numbers for the a-summation 
in ( |l.2[) . As explained in this problem can only be solved if the notion of "negative- 
energy states" is given up and replaced by a causality principle for the Dirac sea; this 
gives a unique definition of P in terms of a power series in the external potential. Our 
task is to convert this formal definition into explicit formulas for the fermionic projector 
in position space. The basic technique is to construct solutions of the inhomogeneous 
Klein-Gordon and Dirac equations and to show that these solutions coincide with the 
contributions to the perturbation expansion of P. For the contribution to first order in the 
external potential, a similar technique was already used in || , which also contains a general 
discussion of the method. In the following, we will first generalize this technique to higher 
order perturbation theory. Then we will explicitly sum up the light-cone expansions of 
all Feynman diagrams, which will finally yield exact formulas for the light-cone expansion 
of the fermionic projector without the restriction for the external potential to be (in any 
sense) "small." We shall use the notation, the definitions, and the results of Q throughout. 
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Since the method of 0] had to be refined considerably for the analysis of higher order 
Feynman diagrams, we will develop the light-cone expansion from the very beginning. 
Thus this paper can be considered as being independent of [||. Nevertheless, the more 
elementary approach in Q is a preparation and might be helpful for the understanding. 
We will use the so-called residual argument to deduce the light-cone expansion of the Dirac 
sea from that of the advanced and retarded Green's functions. This allows us to bypass the 
explicit Fourier transformations in [0]. However, the residual argument has its limitations; 
making it mathematically precise leads to the decomposition of the Dirac sea into a causal 
and a non-causal contribution. 

In the remainder of this section, we specify our problem in mathematical terms. Since 
a realistic physical system consists of several types of fermions, we describe empty space 
by the fermionic projector of the vacuum, which was introduced in |Q] as the direct sum of 
/ > 1 Dirac seas. The chiral asymmetry matrix X and the mass matrix Y are considered 
as a-priori given. The reader who is only interested in the light-cone expansion of a single 
Dirac sea may specialize to / = 1 and X = 1 = Y. On the wave functions, we consider 
the indefinite scalar product 



<^|$> = /_ ^ *z0*0 drx (1.4) 
i=i 



IR 4 



with the adjoint spinor ^ = ^*7°. Similar to the interaction is described by 

a perturbation B of the Dirac operator. We allow B to be composed of chiral and 
scalar /pseudoscalar potentials, 

B(x) = X l4r(x) + Xr4l( x ) + + n 5 ^(x) , (1.5) 

where Xl/r = |(1 -F 7 5 ) are the chiral projectors and where we use a matrix notation in 
the Dirac sea index, 

A L/R = { A L/Rm)l,m=l,...,f , $ = ($ m )/,m=l,...,/ , 3 = (H m )i, m =l,...,/ (1.6) 

(so the potentials may be non-diagonal on the Dirac seas and thus describe a general 
interaction of all the fermions). Furthermore, the perturbation B shall be Hermitian 
with respect to the scalar product ( |1.4| ). We assume the Dirac operator to be causality 
compatible with X, i.e. 

X*(i0 + B-mY) = (i(f) + B-mY)X . 

This assumption is crucial; if it was violated, unbounded line integrals would occur in 
the light-cone expansion, making it impossible to carry out the sum over all Feynman 
diagrams (see the calculations in |J for more details). The form of the chiral decomposition 
Xl4-r + Xr4l ls use f m because, as we will see later, the left and right handed components 
of the fermions couple to Al and Ar, respectively. An interesting feature of our system is 
that, as a consequence of the non-diagonal form (1.6) of the potentials on the Dirac seas, 



the potentials and the mass matrix do in general not commute with each other, 

[A L (x),A L/R (y)]^0 , [A L (x)My)]¥=0 , [A L (x),Y] ? , . . . . (1.7) 

Compared to the situation in Q, the bilinear potential Hj^a^ is missing in (|0|). It 
leads to complications when the sum over all Feynman diagrams is carried out. These 
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complications are not serious, but in order to keep the expansion reasonably simple, the 
bilinear potential was left out. Furthermore, we do not consider the gravitational field. 
The reason is that the higher order contributions in the gravitational potential become 
more and more singular on the light cone. This leads to technical problems which we will 
not deal here. Despite these simplifications, the considered ansatz for B includes arbitrary 
left and right handed Yang-Mills potentials and is general enough for a description of 
e.g. the interactions of the standard model. The fermionic projector in the presence of 
external fields, P(x, y), is defined via the perturbation series in Q, which is a formal sum 
of operator products of the form 

oo 

P = EE const(n, . . .) C n ,... B C n ^„, B ■■■ B C ,... , 

71=0 ••• 

where the factors Cj ... coincide either with the spectral projectors k, p, or with the Green's 
function s (the index '. . . ' is a short notation for the different configurations of these 
factors). 

In the language of Feynman diagrams, the perturbation series for P only consists of 
tree diagrams. These tree diagrams are all finite; this is not difficult to prove if we assume 
suitable regularity of the potential: 

Lemma 1.1 Let (Cj), < j < n, be a choice of operators Cj = k, p, or s. If the external 
potential ( \1.5] ) is smooth and decays so fast at infinity that the functions B(x), x % B(x), 
and x % x°B[x) are integrable, then the operator product 

(C n BC n ^B---BC )(x,y) (1.8) 

is a well-defined tempered distribution on IR 4 x IR 4 . 

Proof: Calculating the Fourier transform of (|L^) gives the formal expression 

\ f d 4 p! f d A p n „i 

xC n (q 2 ) B(qi -Pn-l) C n -i(j> n -l) B(Pn-l ~ Pn-2) ■ ■ ■ Ci(pi) B(pi - qi) C (qi) , (1.9) 

where we consider the Cj as multiplication operators in momentum space and where B 
denotes the Fourier transform of B. It is more convenient to work in momentum space 
because the operators Cj then have a simpler form. We will show that M(q2,qi) is a 
well-defined tempered distribution; our Lemma then immediately follows by transforming 
back to position space. 

The assumptions on B yield that B is C 2 and has rapid decay at infinity, i.e. 

sup \q h ■ ■ ■ q in d K B(q)\ < oo 
<?elR 4 , |k|<2 

for all n, tensor indices h, . . . ,i n , and multi-indices k (so k = (k , . . . , k p ), \k\ := p). As 
is verified explicitly in momentum space, the distributions k, p, and s are bounded in 
the Schwartz norms of the test functions involving derivatives of only first order, more 
precisely 

|C(/)| < const ||/||4,i with C = k,p, or s and / G S, 
where the Schwartz norms are as usual defined by 

P,q = , m & X SU P W 9jf{x)\ . 

\I\<p, \J\<1 xeiR 4 
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As a consequence, we can apply k, p, and s even to functions with rapid decay which 
are only C 1 . Furthermore, we can form the convolution of such functions with k, p, or 
s; this gives continuous functions (which will no longer have rapid decay, however). A 
convolution decreases the order of differentiability of the functions by one. 
We consider the combination of multiplication and convolution 

F{P2) := J-^f(P2-Pi)C( Pl )g(p 1 ) , (1.10) 

where we assume that / € C 2 has rapid decay and g G C 1 is bounded together with its 
hrst derivatives, ||<?||o,i < oo. For any fixed P2, the integral in ( [1.10 ) is well-defined and 
finite because f(p2 — • ) g{) is C 1 and has rapid decay. The resulting function F is C 1 and 
bounded together with its first derivatives, more precisely 

||F|| ,i < const ||/|| 4 ,2 N|o,i • (1-11) 

After these preparations, we can estimate the integrals in (|0]) from the right to the 
left: We choose two test functions f,g € 5(IR 4 ,(D 4 ^) and introduce the functions 



^i(pi) = / 7^s(pi-3i)Co(gi)0(gi) (1.12) 



A 2 

FjiPj) = J -t|^ B( Pj -pj-i) Cf-ifo-i) ^-ife-i) , Kj<n. (1.13) 

The integral ( 1.12| ) is of the form Ql.iop and satisfies the above considered assumptions on 
the integrand. Using the bound ( |1.11| ), we can proceed inductively in ( |1.13|) . Finally, we 
perform the (^-integration, 



M(f, 9 ) = J ^f(q2)C n (q 2 )F n ( q[ 



We conclude that M is a linear functional on <S(IR 4 ,(D 4 ^) x <S(IR 4 ,© 4 ^), which is bounded 
in the Schwartz norm 1 1 . 1 1 4,1 of the test functions. ■ 

Clearly, the existence of the perturbation expansion to every order does not imply the 
convergence of the perturbation series, and we will not study this problem here. Our 
method is to first perform the light-cone expansion of the individual Feynman diagrams. 
For the resulting formulas, it will then be possible to sum up the perturbation series. Since 
the Feynman diagrams are only defined as distributions, we must generalize (1.3) in a way 
which allows for the possibility that A(x, y) is singular on the light cone. 

Def. 1.2 A tempered distribution A(x,y) is of the order 0((y — x) 2p ), p£ Z, if the product 

(y-x)~ 2p A(x,y) 

is a regular distribution (i.e. a locally integrable function). It has the light-cone expan- 



sion 



A(x,y) = y) (1.14) 



3=9 
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with g £ Z if the distributions A^(x,y) are of the order 0((y — x) 2j ) and if A is approxi- 
mated by the partial sums in the way that 



A(x,y) — A"* (x, y) is of the order 



0((y 



(1.15) 



3=9 



for all p > g. 



The lowest summand A^ 9 \x, y) gives the leading order of A(x, y) on the light cone. If A is 
singular on the light cone, g will be negative. The light-cone expansion (|ll^) of a smooth 
function is recovered as a special case by setting g = and A^(x, y) = (y — x) 2j Aj(x, y). 
Notice that the definition of the light-cone expansion does not include the convergence of 
the infinite sum in ( 1.14j ), which is only a convenient notation for the approximation of 
A(x,y) by the partial sums (1.15). 



2 The Light-Cone Expansion of the Green's Functions 

In this section, we shall perform the light-cone expansion for the advanced and retarded 
Green's functions s v , s A , which are defined in [ffl] by the perturbation series 

00 oo 
k=0 k=0 

These perturbation expansions are causal in the sense that s v (x,y) and s A (x,y) only 



depend on the potential B in the "diamond" (L^ nLj)U (L A n L v ), where 



Ll = {y\(y-x) 2 >0, y°-x°>0} , L A = [y \ (y - xf > 0, y° - x° < o} (2.2) 

denote the future and past light cones around x, respectively. Since this "diamond" is 
(for fixed x and y) a bounded region of space-time, the assumptions of Lemma 1.1 on the 
decay of the external potential at infinity are not necessary in this section; it suffices to 
assume that B is smooth. Furthermore, the chiral asymmetry matrix X and the causality 
compatibility condition will not be needed in this section. 

In order to get a first idea of how the light-cone expansion can be carried out, we look 
at the free advanced Green's function of a single Dirac sea in position space: It is 
convenient to pull the Dirac matrices out of by setting 

sL{x,y) = (i$ x + m) S^ 2 (x,y) , (2.3) 

where 5 V 2 is the advanced Green's function of the Klein-Gordon operator, 



SUx,y) = M.„ / ^ -^—23—0 e-*C»-v) . (2.4) 



d A p 1_ 

o<e^»o J (27r) 4 p 2 -m 2 - iep° 

This Fourier integral can be calculated explicitly; we expand the resulting Bessel function 
J i in a power series, 



47r y/m^ (y — xy 
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= ~ S((y- i f) e(y° -x°) 

(0 is the Heaviside function O(t) = 1 for i > and Q(t) = otherwise). This computation 
shows that S^ 2 (x,y) has a <5((y — x) 2 )-like singularity on the light cone. Furthermore, one 
sees that is a power series in m 2 . The important point for the following is that the 
higher order contributions in m 2 contain more factors (y — x) 2 and are thus of higher order 
on the light cone. More precisely, 



d 



dm 2 J m \ m 



S^ 2 , 2=0 (x,y) is of the order 0((y - x) 2n - 2 ). (2.6) 



According to ( |2.3| ), the Dirac Green's function is obtained by computing the first partial 
derivatives of ( |2.5| ). Therefore s^ n {x,y) has a singularity on the light cone which is even 
~ 5'((y — x) 2 ). The higher order contributions in m are again of increasing order on the 
light cone. This means that we can view the Taylor expansion of (|2.3|) in m, 



1 / d \ n 

71=0 ' 

as a light-cone expansion of the free Green's function. Our idea is to generalize this formula 
to the case with interaction. More precisely, we want to express the perturbed Green's 
function in the form 



/ d xn 

r n {x,y) 

n=0 



S v (s,y) = $>„(*, y) \ ^ 2|m2=0 (x,y) (2.7) 



with factors F n which depend on the external potential. We will see that this method is 
very convenient; especially, we can in this way avoid working with the rather complicated 



explicit formula (|2.5|), Apart from giving a motivation for the desired form (2/7) of the 
formulas of the light-cone expansion, the mass expansion ( |2.5| ) leads to the conjecture that 
even the higher order contributions in the mass to the perturbed Green's functions might 
be of higher order on the light cone. If this conjecture was true, it would be a good idea 
to expand the perturbation expansion for s with respect to the parameter m. Therefore 



our basic strategy is to first expand (2.1) with respect to the mass and to try to express 
the contributions to the resulting expansion in a form similar to Q2.7| ). 

The expansion of (|2.l| ) with respect to m gives a double sum over the orders in the 
mass parameter and in the external potential. It is convenient to combine these two 
expansions in a single perturbation series. For this, we first write the mass matrix and 
the scalar /pseudoscalar potentials together by setting 

Y L (x) = Y - — (*(x) + *3(x)) , Y R (x) = Y - - - »H(x)) . (2.8) 

m m 

The matrices Y L / R (x) are called dynamic mass matrices; notice that Y£ = Yr. With this 
notation, we can rewrite the Dirac operator in the form 

i^ + B-mY = i@+B with (2.9) 
B = XL(4 R -m Y R ) + X R (4l ~mY L ) . (2.10) 
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For the light-cone expansion of the Green's functions, we will always view B as the per- 
turbation of the Dirac operator. This has the advantage that the free theory consists only 
of zero-mass fermions; the Green's functions of the free Dirac operator have the simple 
form 

s v (x,y) = i$> x S^ 2=0 (x,y) , s*(x,y) = i<f) x S£ 2=0 (s,y) (2.11) 

(to be very precise, we should write s y = I^Sq (8)1, where 1 = (l m )i,m=l,...J is the identity 
on the Dirac seas). The Green's functions with interaction are given by the perturbation 
series 

oo oo 
~V = £ ( _ s V B) fc s V ^ ~A = £(_ s A£)fc s A _ (212) 

k=0 k=0 

The constructions of the following subsections are exactly the same for the advanced 
and retarded Green's functions. In order to treat both cases at once, we will in the 
remainder of this section omit all superscripts ' v ', ' A '. The formulas for the advanced and 
retarded Green's functions are obtained by either adding ' v ' or ' A ' to all factors s, S. 



2.1 Inductive Light-Cone Expansion of All Feynman Diagrams 

In this subsection, we explain how the individual contributions to the perturbation expan- 
sion ( 2.1 2| ) can be written similar to the right side of fl2.7D as a sum of terms of increasing 
order on the light cone. For the mass expansion of S m 2, we set a = m 2 and use the 
notation 

S {1) = (£)^ a | 0=0 . (2.13) 

In preparation, we derive some computation rules for the S^ 1 ': S a satisfies the defining 
equation of a Klein-Gordon Green's function 

(-□x - a) S a (x,y) = 5 4 (x-y) . 

Differentiating with respect to a yields 

-D x S®(x,y) = 5 lfl 5\x-y) + I S^ l \x,y) , I > 0. (2.14) 

(For I = 0, this formula does not seem to make sense because S^ 1 ^ is undefined. The 
expression is meaningful, however, if one keeps in mind that in this case the factor I is 
zero, and thus the whole second summand vanishes. We will also use this convention in 
the following calculations.) Next, we differentiate the formulas for S a in momentum space, 

SaiP) = "5 ^—^7 , S£(p) = -= l — -= , (2.15) 

p 2 - a- iep° p 2 -a + iep° 

with respect to both p and a. Comparing the results gives the relation 

^S a ( P ) = -2 Pk ±S a {p) , 

or, after expanding in the parameter a, 

-§^S^(p) = -2p k S^ l+l \p) , l>0. (2.16) 
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This formula also determines the derivatives of in position space; namely 



dx k K ' J (2vr) 4 

(2.16) * r d 4 p d -ip( x -y) _ _£ / ^ c(J-i)(„l JL P -^- 

2j(2TT^dp k [P> 2 J (2vr) 4 W 



•y) 



= (*,!/) , Z>1. (2.17) 

We iterate this relation to calculate the Laplacian, 



2S ( - [ - 1 \x,y) - ±(y- x fs( l - 2 \x,y) 



I > 2. 



After comparing with ( |2.14 ), we conclude that 

(y-x) 2 S (l) {x,y) = -41 S (l+1) {x,y) , l>0. (2.18) 

Finally, S^ l \x,y) is only a function of (y — x), which implies that 

d -SV\x,y) = Os®(x,y) , I > 0. (2.19) 



dx k ' dy k 

The following lemma gives the light-cone expansion of an operator product which is 
linear in the external potential. We will later use it for the iterative light-cone expansion 
of more complicated operator products; in this case, the potential will be a composite 
expression in B and its partial derivatives. In order to avoid confusion then, we denote 
the external potential by V . 

Lemma 2.1 (light-cone expansion to first order) The operator product V 
with l,r > has the light-cone expansion 

(S® V S<ri)(x,y) 

£ / a l (l-aY(a-a 2 r(n n V) lay+{1 . a)x daS( n+l+r+1 Hx,y) . (2.20) 



Proof: The method of the proof is to first compute the Laplacian of both sides of ( 2.20Q . 
The resulting formulas will be of similar structure, which allows us to proceed inductively. 



On the left side of ( |2.20| ), we calculate the Laplacian with the help of ( 2.14 ) to 

-n x (S®VSV)(x,y) = 8 lt0 V(x) S<r\x,y) + I (S^ V S^)(x,y) . (2.21) 

The Laplacian of the integral on the right side of ( |2.20| ) can be computed with ( |2,17j ) 
and ( gig ), 

-U x f a 1 (1 - af (a - a 2 ) n (U n V)\ ayH ^ a)x da S^ n+l+r+1 \x, y) (2.22) 
a 1 (1 - a) r+2 (a - a 2 ) n {U n+l V)\ ay+(l „ a)x da S^ l+r+1 \x,y) 
a 1 (1 - a) r+1 (a - a 2 ) n {d k U n V)\ ay+{l . a)x da (y - x) k S^ n+l+r \x,y) 
+(n + / + r + l) f\ l (I - a) r {a - a 2 ) n {n n V)\ ay+{l _ a)x da S^+^Xx^y) . 

J 
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In the second summand, we rewrite the partial derivative as a derivative with respect to 
a and integrate by parts, 



a 1 (1 - a) r+1 (a - a 2 ) n (d k U n V)\ ayHl . a)x da (y - x) k 
= j a 1 (1 - a) r+l (a - a 2 ) n — {U n V)\ ay+(l . a)x da 

-S n ,o Si,o v ( x ) - ( n + / C 1 - «) r+2 (« " « 2 ) n_1 (a n V0|ay+(i-a)* da 

J 

+(n + r + 1) jf a' (1 - a) r (a - a 2 ) n (□ n F)| a2/+(1 _ a)a; da 
-n [ a 1 {I- a) r+2 (a - a 2 )™" 1 (□ n F) K+(1 _ a):r da 

« 

+ (n + / + r + 1) / a' (1 - a) T (a - a 2 ) n (n n V) ]ay+{1 „ a)x da 

J 

-I f a 1 ' 1 (1 - a) r (a - a 2 )" (□ ri V) K+( i_ a):E da . 

J 

We substitute back into the original equation and obtain 

+/ /' a'" 1 (1 - a) r (a - a 2 ) n {U n V)\ ay+(1 _ a)x da S^ +l+r \x,y) 

j 

- f a 1 (1 - af +2 (a - a 2 ) n (U n+l V\ ay+{l _ a)x da S^ n+l+r+1 \x,y) 
x l (1 - a) r+2 (a - a 2 )"" 1 (n n V)\ ay+(1 „ a)x da S^+ l+r \x,y) . 



-n I a" 



After dividing by n! and summation over n, the last two summands are telescopic and 
cancel each other. Thus one gets 

-□El (1 - «)' (« - « 2 )" (□'^Wa-a)* A* S^ n+l+r+1 Hx,y) 

= S lj0 V(x)S^(x,y) 

oo -, i 

+^E-/ ^-\l-aY(a-a 2 r(n n V) lay+{1 _ a)x daS^ +l+r Hx,y) . (2.23) 

n=0 n - J ° 

We now compare the formulas ( |2.21| ) and ( 2.23j ) for the Laplacian of both sides of 
( |2.20| ). In the special case I = 0, these formulas coincide, and we can use a uniqueness 
argument for the solutions of the wave equation to prove ( 2.2C| ): We assume that we 
consider the advanced Green's function (for the retarded Green's function, the argument 
is analogous). For given y, we denote the difference of both sides of ( |2.20|) by F(x). Since 
the support of F(x) is in the past light cone x E Ly, F vanishes in a neighborhood of the 
hypersurface Tt = {z £ IR 4 | z° = y° + 1}. Moreover, the Laplacian of F is identically zero 
according to ( |2.21| ) and ( |2.23| ). We conclude that 

UF = and F\ n = d k F\ n = . 
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Since the wave equation has a unique solution for given initial data on the Cauchy surface 
H, F vanishes identically. 

The general case follows by induction in I: Suppose that (|2.20 ) holds for given I (and 



arbitrary r). Then, according to ( 2.21| ), ( |2.23 ), and the induction hypothesis, the Laplacian 



of both sides of fl2.2Cp coincides for I = 1 + 1. The above uniqueness argument for the 
solutions of the wave equation again gives ( |2.20[ ) . ■ 



We recall for clarity that, according to (|2.6|), the higher a-derivatives of S a (x,y) are of 
higher order on the light cone. Thus the summands in ( 2.2C| ) are of increasing order on 



the light cone, and the infinite sum makes mathematical sense in terms of Def. |1.2| via the 
approximation by the partial sums ( |1.15| ). 



The representation ( p. 20 ) of an operator product as an infinite series of line integrals 



has some similarity with the formal light-cone expansion g, Theorem 3.3]. In order to 
make it easier for the reader to see the connection between these two expansions, we 
briefly discuss the analogy and the differences: First of all, we are here considering the 
Green's functions instead of the negative-energy solutions of the Klein-Gordon equation. 
The causality of the Green's functions (i.e. supp S(x, .) C L x ) simplifies the construction 
considerably. We could use it to avoid the explicit Fourier transformations of the proofs 
in Q; furthermore, it makes the resummation of non-local contributions unnecessary. In 
this paper, the complications related to the non-causality of the negative-energy states 
will reappear in the light-cone expansion of the Dirac sea in Section ||[ We point out 
that the light-cone expansion in [g| is more general in the way that the mass parameter 
a = m 2 need not be zero. This is sometimes more convenient, because then the mass is 
just a parameter of the Green's functions instead of occurring in the line integrals of the 
light-cone expansion. With our concept of dynamic mass matrices, however, an expansion 



around a = is more appropriate. The important generalization in ( 2.20 ) is that the two 
factors S on the left side may be derivatives with respect to a. On the right side of equation 
( 2.20D , this is taken into account by the additional factors a 1 (I — a) r in the integrand and 



by the higher order (n + l + r + 1) of the a-derivative of S. This generalization is the basis 
for the following iterations. 

Lemma can be used for the light-cone expansion of more complicated operator 
products. To explain the method, we look at the simplest example of three factors S^ ' 
and two potentials V, W, 

( 5 (0) y 5 (Q) w sW)(x,y) = fd 4 zS ( -°\x,z)V(z)(S^WS^)(z,y) . (2.24) 



Having split up the operator product in this form, we can apply Lemma £j] to the factor 

= E ^/ d * z S {0 \x,z) fa) j\a - a 2 ) n (D n W) lay+{1 _ a)z da} S^ +1 \z,y) . 

Now we rewrite the z-integral as the operator product (S^ g y S^°')(x,y), where g y (z) is 
the function in the curly brackets. The y-dependence of g y causes no problems because 
we can view y as a fixed parameter throughout the expansion. Thus we can simply apply 



Lemma |2.1| once again and obtain 

= E ZTZrJ (1 - (3T+ 1 05 - /3 2 ) w / da(a-a 2 r 

„ fib. lb. J(\ JQ 



n 



x n?(v(z) (n n w) 



\ay+(l— a)z 



\z=(3y+(l-p)x 



S^ m+n+2 \x,y) 



The Laplacian could be further computed with the Leibniz rule. Notice that the 
manipulations of the infinite sums are not problematic because the number of terms to 
every order on the light cone is actually finite (the situation would be more difficult if 
we studied the convergence of the sum (1.14), but, as pointed out earlier, the light-cone 
expansion is defined merely via the partial sums). 

We want to iteratively perform the light-cone expansion of the operator products in 
( |2.12| ). This is not directly possible with the method just described, because (2.12) contains 
the Dirac Green's function s (instead of S). We must think about how to deal with this 
problem. Relation ( 2.1 1] ) allows us to replace the factors s by S, but this gives additional 
partial derivatives in the operator product. These derivatives can be carried out after each 
iteration step using the Leibniz rule and the differentiation rule ( 2.17| ). In the simplest 
example, we have 

(s^VS^)(x,y) = (i@ x )(S^VSW)(x,y) 



E 



n=0 



n\ Jo 



[a 



,2\n 



a 



(n n V) lay+{1 _ a)x S^ +1 \x,y) 



; E 

n=0 



1 

n! Jo 



(1-a) (a - a 2 ) n {$n n V) 



\ay+(l— a)x 



S {n+1) (x,y) 



Z n=0 n - J0 

The only problem with this method is that the partial derivatives might hit a factor S^°\ 
in which case rule ( |2.17| ) cannot be applied. In order to resolve this problem, we extend 
our constructions in a way which allows us to use all previous formulas also in this special 
case. For this, we take ( [2.17 ) as the defining equation for (y — x)k S^ l \x,y), 



(y-x) k S < - l \x,y) 



d 
dx k 



(2.25) 



(notice that S^ 1 ^ itself remains undefined, only the combination (y — x)k S^ -1 )^, y) 
makes mathematical sense as the partial derivative of the distribution S^). With this 
definition, the computation rule ( 2.18 ) becomes also valid for I = — 1: According to ( |2.18| ), 
the distribution (y — x) 2 S^°'{x,y) vanishes identically, and thus 



= n x ((y -x) 2 S^(x,y) 

= (D x (y-x) 2 )s^(x,y) - 4(y- X y : JL-S<-°\x,y) + (y - x) 2 n x S^ (x,y) . 

In the last summand, we substitute ( |2.14j) , which gives —{y — x) 2 6 4 (x — y) = 0. Using the 
identity U x {y - x) 



8 and the definition ( |2.25 ), we conclude that 
(y-xfS^ix^) = 4S®(x,y) 



(2.26) 



The following lemma extends the result of Lemma 2A to the case when the right factor 
of the operator product is S^ 1 ^. 
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Lemma 2.2 (light-cone expansion to first order for r = — 1) The operator product 
(5(0 . St" 1 )), / > 0, has the light -cone expansion 

d 4 zSV(x,z)V(z)(y-z) k S(- 1 \z,y) = £ - / a 1 (1 - a)" 1 (a - a 2 ) n 
J ^ n\ Jo 

xn n AV(z)(y-z) k ) lz=ay+{1 _ a)x daSW(x,y) . (2.27) 

Proof: We deduce the light-cone expansion for the left side of Q2.27 ) from (2.20) by 
pulling one y-derivative outside, 



d*zS«\x,z)V(z)(y-z) k S^(z,y) (2 ' 25) =» -2 y S ^y) 



(2.19), (2.17) 



-2 p E^/ «' (« - « 2 ) n da 5(™ +/+1 )(x, y) 

° n=0 ' ^ 
oo -, ! 

^ / a ' ( a _ a 2)« (□"XOlcnH-d-a)* <*a (y - x) k S^ n+l \x,y) (2.28) 

rc=0 ,/0 

oo -. „! 

" 2 E-/ « m («-« 2 f (^□^)| Q , +(1 _ a)a;( ia^ + ' +1 )(x,y) . (2.29) 



In ( 2.29 ), we substitute the identity 

d k n n v(z) - 1 



U n z +l (V(z) (y-z) k ) + 



1 



D n z +1 V(z)) (y-z) k . (2.30) 



2(n + l) z v w v ' ' 2(n + l) 

The contribution of the second summand in fl2.30| ) cancels all the summands n = 1,2, 
of (|2.28| ), and we get 

d 4 z S® (x, z) V(z) (y - z) k (z, y) = [ a 1 V ]ay+{1 _ a)x da (y - x) k S® (s, y) 



oo -. r i 

+ Et-tw/ ^ 1 ( a -^ra^' 1 (^W(»-^)*)|^(i^)-^^ 1+,) ( a! .i') • 

n=0 ^ ^ J0 

After shifting the summation index, this coincides with ( 2.27] ). ■ 

Notice that the pole of the factor (1 — a) -1 for a = 1 does not lead to problems in (2.27): 
In the case n = 0, it disappears since (1 — a) -1 (y — z) = y — x; for n > 0, it is compensated 
by the zero of the factor (a — a 2 ) n . 

After these preparations, we come to the light-cone expansion of general Feynman 
diagrams. For the line integrals of Lemma 2.1 and Lemma 2.2, we introduce the short 
notation 
rv 

(2.31) 



[l,r \n]dz f(z) := [ da a 1 (1 - a) r (a - a 2 ) n f(ay + (1 - a)x) 
Jo 

If I, r, and n all vanish, we sometimes omit the bracket, i.e. 



dz f(z) := / da f(ay + (1 — a)x) 



Furthermore, we abbreviate the following products with multi-indices 



dl :-- 



d d 



(y - x) J := (y - x) n ■ ■ ■ (y - x) n , -y J := -y n ■ ■ ■ ^ , 



dz^ dzii 
where J = (j'i , • • • 
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Theorem 2.3 (inductive light-cone expansion) The light-cone expansion of the k th 
order contribution {(—a B) k s)(x,y) to the perturbation series /\2.1i ) can be written as an 
infinite sum of expressions, each of which has the form 

Xc C(y- x) 1 [ y [h, ri | m] dz x dil Ul\ V% ( Zl ) [ y [l 2 ,r 2 I n 2 ] dz 2 d{\ U% V% C2 (z 2 ) 

•••f [lk,r k \n k }dz k diln p z tV^ Ck (z k )r J S( h Hx,y) . (2.32) 

J Zu—l 



In this formula, C is a complex number and the parameters l a , r a , n a , and p a are non- 
negative integers; the indices c and c a can take the two values L or R. The functions 
Vj c coincide with any of the individual potentials in (2.1C) with chirality c a , i.e. 

y(a) _ ^ ^ n yjfacfo case \J a \ = l) or 
Vjf> = mY Ca (in which case \J a \ = 0) . (2.33) 
The chirality c a of the potentials is determined by the following rules: 

(i) The chirality c\ of the first potential coincides with the chirality of the projector \c- 

(ii) The chirality of the potentials is reversed at every mass matrix, i.e. 

, f coincide if Vc„ = A Cn 



The tensor indices of the multi-indices in (2.3i) are all contracted with each other, ac- 
cording to the following rules: 

(a) No two tensor indices of the same multi-index are contracted with each other. 

(b) The tensor indices of the factor j J are all contracted with different multi-indices, in 

the order of their appearance in the product ( 2.32[ ) (i.e., for J = (ji,...,j/) and 
1 < a < b < I, the multi-index with which j a is contracted must stand to the left of 
the multi-index corresponding to ji). 

The parameter h is given by 

k 

2h = k-l-\I\ + + 2 Pa) ■ ( 2 - 34 ) 



0=1 



The number of factors (y — x) is bounded by 

k 

|/| < fc + l-£|/o| • (2.35) 

a=l 

Basically, this theorem states that the light-cone expansion of the A: th order Feynman 
diagrams can be written with k nested line integrals. Notice that, according to ( |1.7| ), 
the potentials V^ a \z a ) do in general not commute with each other, so that the order 
of multiplication is important in ( 2.32; ). In order to avoid misunderstandings, we point 
out that the derivatives d\ a a and do not only act on V^ a \z a ), but also on all the 
following factors V^ a+l \z a+ \), V^ a+2 \z a+2 ),. . . (note that the variables z a +i, z a+2) ... 
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implicitly depend on z a via the inductive definition of the line integrals). Clearly, these 
derivatives could be carried out further with the Leibniz rule, but it is easier not to do 
this at the moment. The restrictions (a) and (b) on the possible contractions of the 
tensor indices were imposed in order to avoid an abuse of our multi-index notation. More 



precisely, (a) prevents factors (y 



in (y 



an unnecessary large number of 7- 



matrices in j J , and "hidden" Laplacians in the partial derivatives d\ a a . The rule (b), on 
the other hand, prevents factors (y — x) 2 and hidden Laplacians in combinations of the 
form (y — x)i (y — x)j 7* 7 J and dijVj 7* 7 J , respectively. Our ordering condition for 
the 7-matrices is just a matter of convenience. Relation ( 2.34[ ) is very useful because it 
immediately tells for any configuration of the line integrals and potentials in Q2.32 ) what 
the corresponding order on the light cone is. 
inequality 



Notice that ( gjg ) and ( ggg ) imply the 



h > -l + J2\Ia\+Pa 



(2.36) 



a=l 



Especially, one sees that h > — 1. In the case h = -1, (fl34|) yields that \I\ > 0, so that 
(2.32) must contain at least one factor (y — x). Therefore the factor S^ h ' in ( p. 32 ) is always 



well-defined by either (|2.13[) or (|2.25|) . 

We point out that, although the total number of summands ( 2.32j ) is infinite, the 
number of summands for any given value of the parameter h is finite. This is clear because, 
for fixed h, the relations ( |2.34| ) and ( 2.35 ) only allow a finite number of possibilities to 
choose the parameters |/|, \I a \, and p a ; thus one only gets a finite combinatorics for all 
expressions of the form ( ^32|) . Since, according to ( |2.6| ), the contributions for higher values 
of h are of higher order on the light cone, we conclude that the number of summands ( 2.32| ) 
is finite to every order on the light cone. Therefore the light-cone expansion of Theorem 



makes mathematical sense in terms of Def. 1.2 



Proof of Theorem \2. We proceed inductively in k. For k = 0, the assumption is true 
because we can write the free Dirac Green's function in the required form ( |2.32 ), 



s(x,y) 



(2.11 ),(|2.25| ) 



(XL + XR) \ (y - xVlj S { 1} (x,y) 



(2.37) 



The conditions (i), (ii), (a), (b), and the relations fl2.34|) , ( 2.35| ) are clearly satisfied. 
Assume that the theorem holds for a given k. With the formula 



((sB) k+1 s)(x,y) 



(2.11) 



otz S(°\x,z)B(z) ((-a B) k s)(z,y) 



(2.38) 



we can express the (k + l) st order contribution to the perturbation series ( |2.12| ) in terms 
of the fe th order contribution. We must show that ( |2.38 ) can again be written as a sum of 
expressions of the form ( |2.32| ) (with k replaced by k + 1), and that (i), (ii), (a), (b), and 
( |2.34| ), ( 2 .35) ) are satisfied. This is done in several construction steps: 



1) Chiral decomposition: We substitute the induction hypothesis ( 2.32j ) into ( 2.38] ). This 
gives a sum of expressions of the form 

Ci^J d'z S(°\x,z) ((y - z) 1 B{z) Xc j\h,n \ ni] dz x d{\ n pi vj^z,) 

(2.39) 



v 

Zk-l 



[i k , rk 1 n k ] dz k aii ° Pfc v J k 2cS z *) ~f J S(h) ^y) ■ 
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We insert the special form of the potential B, ( [2.10 ), and expand. Using the commu- 
tation rule Yxl/r = XRjLl' \ we bring all chiral projectors to the very left, where they 
can be combined with the formula XcXd = ^cd Xc to a single chiral projector. Next, 
we bring the 7-matrices of B to the right and write them together with the factor / y J 
in ( 2.39| ) (notice that the Dirac matrices commute with the potentials Vc " , which 
are only non-diagonal in the Dirac sea index). Denoting the individual potentials of 
the factor B in (|2.39| ) by Vj®^, we thus get for ( p. 39 ) a sum of expressions of the 
form 



XcCi^J d i zS(°\x,z)l l (y-z) I V^ co (z) Jjl 1 ,r 1 \n 1 ]dz 1 d£n» 

[i k ,n I nk ] dz k dii n Pk v^ Ck {z k ) A s^(z, y ) . (2.40) 



y 

Zk-l 



The chiral decomposition (2.10) and the induction hypothesis (i) yield that the 



chirality cq coincides with c, whereas c\ coincides with cq if and only if V c 



(0) 



CO 



Since the chiralities C2, C3,. . . satisfy the induction hypothesis (ii), we conclude that 
the rules (i) and (ii) are also satisfied in ( |2.40| ) (after relabeling the indices in an 
obvious way). The chirality of the potentials will not be affected in all the following 
construction steps; to simplify the notation, we will omit the indices c a from now 



on. 



2) Light-cone expansion with Lemma 2.1, Lemma 2.i: Since y can be considered as a fixed 



parameter, we can in ( 2.40 ) apply either Lemma 2A or Lemma 2^2 with V given by 
the expression in the curly brackets, 

00 1 ry 

(2.40) = XcCi$ x Y,-\ / %h\n]dz 

n=0 n - Jx 

]J ((y - z) 1 < } W fjh,n I m] d Zl flg v£\ Zl ) 



x □ 



V 

Zk-l 



[k,r k I n k ] dz k dil n Pk V£\z k ) ) r J S^ n+h+1 \x,y) 



(2.41) 



3) Computation of the Laplacian We carry out the z-derivatives in ( |2.41 ) inductively 
with the Leibniz rule. Each derivative can act either on the factors (y — z) 1 or on 
the functions V^ a \ In the first case, one of the factors (y — z) disappears. Thus we 
get a sum of expressions of the form 

Xc C i$ x f V [0, h\n]dz(y- z) 1 D§° V%\z) f [l u r x \ m] d Zl d% D Pl V^{z x ) 

Jx J z 



y 

I 

Zfc-1 



l k ,r k I n k ] dz k dil V£\z k ) 7 J S (n+/l+1) (x, y) 



Jk 



(2.42) 



with 111 < 1/1 and 



2n = \I\ - \I\ + |J | + 2p . (2.43) 

We can assume that no tensor indices of d^. are contracted with each other (otherwise 
we rewrite the corresponding partial derivatives as additional Laplacians). Then all 



the partial derivatives d z in (2.42) were generated in the case when one derivative 
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of a Laplacian D 2 in ( 2.41 ) hit a factor {y — z) whereas the other derivative acted 
on the V ^ . Thus the number of factors (y — z) which disappeared by carrying out 
the Laplacians in ( |2.41| ) is larger or equal than the number of partial derivatives d z , 

\I\ ~ \I\ > \Io\ • (2.44) 
4) Extraction of the factors (y — x): In ( 2.42; ), we iteratively apply the identity 

(2.45) 



[0, r | n] dz (y — z) ■ ■ ■ = (y — x) / [0, r + 1 1 n] dz ■ ■ ■ 

J X 

This gives (k + 1) nested line integrals of the form 

(2.42) = Xc C i@ x (y - x) 1 S&(x, y) f % \l , r \ n ] dz ^ D Po V^\z ) 



y 



[Z fc> r fc |n fc ] dz k d I z lU Pk V^ ) {z a ) 1 



with 







r = h+\I\ 



< 2h = 2(n + /t + l) 



(2.43) 



uq = n 

2h + 2 + \I\- \I\ + |/ | + 2 Po 



(2.46) 



(2.47) 
(2.48) 



We can arrange that the parameters Iq, ro, and no are all positive: The only param- 
eter which might be negative is tq; in this case, h = —1, \I\ = 0, and thus rg = — 1. 
The induction hypothesis ( p. 34 ) yields that |/| > 0. Thus |/| > \I\, and relation 
( |2.43| ) gives that (n =)n > 0. Therefore we can apply the identity 

[l ,r | no] = [Iq + 1, r + 1 1 n - 1] 

to make all the parameters in this bracket positive. 

5) Computation of the partial derivative $ x : The x-derivative in ( |2.46 ) can act on the fac- 
tors S( h) , (y-x) 1 , or V^{z a ). The first case can be computed with the rules ( |2.17| ) 
or ( |2.25| ); it decreases h by one and gives one additional factor (y — x). In the second 
case, one factor (y — x) disappears, and thus \I\ is decremented. The last case can 
be handled with the rule 



d fy 



dx k 



rv rv 
/ [l,r | n] dz f(z,y) = / [Z,r + l|ra] 

J X J X 



dz k 



f(z,y) 



(2.49) 



which increases |/o| by one. As is immediately verified in each of these cases, equation 
( |2.48| ) transforms into 



2h = 2h + 1 + |/| - |/| + \I \ + 2 Po 
whereas inequality ( 2.44] ) must be weakened to 

|/| < i + |/|_|/ | . 

Finally, we combine the 7-matrix of the factor (j) x with j J . 



(2.50) 



(2.51) 
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After these transformations, the (k + l) st order Feynman diagram consists of a sum of 
terms of the form 

(0), 



Xc C(y- x) 1 / [Z 0) ro | n ] dz dj° n p z ° V}°> (z ) 



v 



[h,n | n k ] dz k dll Ull Vr (z k ) 7 J (x, y) . (2.52) 



z k-l 



Notice that the parameters I a ,p a , a = 1, . . . ,k, were not changed by the above construction 
steps; they are still the same as in the induction hypothesis ( 2.32| ). After renaming the 



indices and the integration variables, ( |2.52| ) is of the required form ( |2.32 ). The conditions 



(a) and (b) for the contractions of the tensor indices, however, will in general be violated. 
Therefore we need two further computation steps: 

6) Simplification of the ^-matrices: If any two of the tensor indices of the factor j J are 
contracted with each other, we reorder the 7-matrices with the anti-commutation 
relations 

{ 7 \ 7^'} = 2g*l (2.53) 

until the corresponding matrices are next to each other. Applying the identity 
7*7i = 41, both Dirac matrices disappear. We iterate this procedure until no tensor 
indices of 7^ are contracted with each other (notice that the iteration comes to an 
end because the number of 7-factors is decreased by two in each step). Again using 
the anti-commutation rule ( 2.53j ), we reorder the Dirac matrices until they are in the 



same order in which the factors to which their tensor indices are contracted appear 
in the product ( 2.52j ). If any two of the 7-matrices are contracted with the same 



multi- index, these 7-matrices are next to each other, and we can use the symmetry 
in the tensor indices to eliminate them both, more precisely 

(y - x)t (y - x)j ■ ■ ■ 7V = (y-x) 2 ---l (2.54) 
%y (a) ---7V = □F (a) ---l . (2.55) 

After all these transformations, condition (b) is satisfied. 

Notice that the parameters \I a \ and p a are in general changed in this construction 
step. More precisely, each transformation ( 2.55| ) modifies the parameters according 
to 

\I a \ -> |4| -2 , Pa -> p a + l . (2.56) 

7) Handling of the new contractions: If any two tensor indices of a factor d l z a a are con- 
tracted with each other, we rewrite the corresponding partial derivatives as a Lapla- 
cian; this changes the parameters |/ a | andp a according to (|2.56 ). If two tensor indices 



of the factor (y — x) are contracted with each other, this gives a factor (y — x) 2 . 
Using the identity ( 2.1S| ) or ( 2.26| ), we inductively absorb the factors (y — x) 2 into 



S^ h \x,y), which transforms h and |/| as 

h -*■ h + 1 , \I\ -»■ \I\ -2 . (2.57) 

After these transformations, condition (a) is also satisfied. 
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After all these construction steps, the (fc+l) th order Feynman diagram is a sum of terms of 
the form ( 2,52j) s atisfying the conditions (a) and (b). It remains to show that the relations 
( gjg) and (|2.35|) remain valid in our inductive construction: As mentioned earlier, the 
parameters I a , p a , a = l,...,k are not changed in the construction steps 1) to 5). In 
the steps 6) and 7), the transformations fl2.56| ) and (2.57) preserve both the induction 
hypothesis (gj^),( gj3g ) and the relations (^50|),(^5l|)7 as is immediately verified. By 



substituting flOpp and ( [2751] ) into (^56|), p37D , we conclude that 



2h 



+ 1)-1-|J| + J2\I a \+2p a 



k 

E 

a=0 



|/| < ( fc +l) + l_^|J a | 



a=0 



Note that this proof is constructive in the sense that it gives a procedure with which 
the light-cone expansion of every Feynman diagram can be carried out explicitly. The 
disadvantage of this procedure, however, is that the resulting formulas become more and 
more complicated to higher order perturbation theory. Therefore it is essential to rearrange 
and collect the contributions of all Feynman diagrams in a way which makes it clear how 
s(x, y) looks like to every order on the light cone. In preparation for this analysis, which 



will be the task of the following Subsection 2/2, we shall now simplify the light-cone 
expansion of Theorem 2.2 a little bit. More precisely, we want to eliminate from ( 2.32j ) 



the partial derivatives in direction of the line integrals, i.e. those derivatives d Zl which 
are contracted with a factor (y — x). We call these derivatives tangential. The following 
combinatorial lemma controls the number of tangential derivatives. 



Lemma 2.4 The contributions ( 2. Si ) to the light-cone expansion of Theorem \2. 3[ satisfy 
for a = 1, . . . , k the inequalities 



l a + n a > t a -l 



and 



r a +n a > ^2t b 

b=a 



(2.58) 



where t a denotes the number of tensor indices of the multi-index I a which are contracted 
with the factor (y — x) 1 . 

Proof: As in the proof of Theorem |2.3| , we proceed inductively in the order k of the 
perturbation theory. For k = 0, the inequalities ( 2.58j ) are trivially satisfied according to 
( 2.37 ). Assume that ( |2.58D is true for a given k. We go through the construction steps 
1) to 7) of Theorem |2.3| and check that the inequalities fl2.58| ) then also hold in (2.52) for 
a = 0, . . . , k. 

We first consider the case a > 0. The parameters l a , r a , and n a remain unchanged in 
all the construction steps of Theorem |2.3| . Furthermore, it is obvious that the parameters 
t a are not affected in the steps 1), 2), 4), and 7). In the steps 3) and 5), the computation 
of the derivatives and $ x might annihilate some of the factors {y — x) which were 
contracted with the factors d^. a ; this may decrease the parameters t a . For the analysis of 
step 6), note that all 7-matrices which are contracted with factors (y — x) stand to the 
left of those 7-matrices which are contracted with the d^, a = 1, . . . , k (this follows from 
the ordering condition (b) in the induction hypothesis and the fact that additional factors 



(y 



7j are only generated during the construction if the partial derivative (j) x hits 



SW i n step 5); in this case, the corresponding 7-matrix stands at the very left in 7^). 
Therefore the commutations of the Dirac matrices do not lead to additional contractions 
between factors (y — x) and d\ a , which implies that the parameters t a remain unchanged 
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in step 6). We conclude that the l a , r a , and n a remain unchanged whereas the t a may only 
decrease, and thus ( 2.5§| ) holds for a = 1, . . . , k throughout all the construction steps. 

It remains to show that the inequalities ( 2.58| ) hold in ( |2,52 ) for a = 0. We first look at 
the situation after step 4) m ( 2.46j ): The values ( |2.47 ) for Iq, ro, and no give in combination 
with ( p.43| ) the equations 



ro + n 



|/| + |/ | + 2p 



1 



h + + \I\ + \M + 2p 



(2.59) 
(2.60) 



Moreover, the number of tangential derivatives to at the first potential is clearly bounded 
by the total number of derivatives there, 



\M > t 



(2.61) 



Furthermore, the total number of tangential derivatives is smaller than the number of 
factors (y — x), 



a=0 

Substituting Q2.44j ) and ( |2.61| ) into ( |2.59j ) yields the inequalities 

lo + n > \I \ +po > t . 



(2.62) 



(2.63) 



In order to get a bound for ro + no, we must distinguish two cases. If h > 0, we substitute 
fl2T44l ) into (j^TBOj) and get with ( |2T62j ) the inequality 



ro + no > \I\ + |/ | + Po > \I\ > ^t a ■ 



(2.64) 



a=0 



In the case h = — 1, ( p. 36 ) shows that \I a \, and consequently also t a , vanish for 1 < a < k. 
Furthermore, ( ggg ) yields that |/| ^ 0. Thus ( pT60| ) and ( ^6T|) , p^62| ) give the bound 

\I\ 1 k 1 1 fc 1 

r + n >/i + — + ~^t a + - to > -^t a + - to , 



a=0 



a=0 



where we used in the last inequality that h + |^|/2 > —1/2 and that all the other terms 
are integers. Since to = ELo*«i we conclude that inequality ( p. 64 ) also holds in the case 
h = -1. 

We finally consider how the bounds ( |2.63j ) and ( |2.64| ) for lo + no and ro + no must 
be modified in the subsequent construction steps. In step 5), the partial derivative X 
may annihilate a factor {y — x), in which case the parameters t a might decrease. On the 
other hand, the partial derivatives ^ may produce an additional factor d ZQ ; in this case, 
ro is incremented according to ( 2.49j ), In step 6), only this additional factor d zo may be 
contracted with (y — x) 1 . Step 7) does not change lo, ro, no, and to. Putting these trans- 
formations together, we conclude that the inequality ( p.63 ) for lo + no must be weakened 
by one, whereas the bound ( |2,64| ) for ro + no remains valid as it is. This gives precisely 
the inequalities ( 2.58|) for a = 0. ■ 
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Theorem 2.5 (partial integration of the tangential derivatives) Every contribu- 



tion (2. Si) to the light cone expansion of Theorem \2.3\ can be written as a finite sum of 
expressions of the form 



Xc C(y-x) K W^(x) [l u r 1 \n 1 }dz 1 W { V{z 1 ) [l 2 , r 2 | n 2 ] dz 2 W^ 2 \z 2 ) 

J X J Z\ 

[l a ,r a | n a ] dz a W (a Hz a ) 7 J S^(x,y) , a < k, (2.65) 



y 

Z a -1 



where the factors are composed of the potentials and their partial derivatives, 

w w = (d K ^ n p^vi ap) ) ■ ■ ■ (d Kh m Pb ^v? 0) 



(2.66) 



with a\ = 1, a/3+1 = bp + 1, bp > ap — 1 (in the case bp = ap — I, is identically one), 

and b a = k. The parameters l a , r a , and n a are non-negative integers, C is a complex 
number, and c = L/R, c a = L/R are chiral indices. The potentials are again given by 
( 2.330 ; their chirality is determined by the rules (i) and (ii) of Theorem 2.3. The tensor 
indices of the multi-indices J, K, J a , and K a are all contracted with each other, according 
to the rules (a),(b) of Theorem iL5 and 

(c) The tensor indices of (y — x) K are all contracted with the tensor indices of the factors 
or 1 J (but n ot with the factors d Ka ). 

We have the relation 

k 

k-l-\K\ + J2(\K a \ + 2p a ) . (2.67) 
a=l 



2h 



Proof: The basic method for the proof is to iteratively eliminate those partial derivatives 
d\l in ([HD which are contracted with a factor (y — x). This is done with the partial 
integration formula 



• f y (2 31) 

(y - x) J / [l,r\ n] dz djf(z) = 

J X 

= <Wn,0 f(y) - ^+n,0 fix) 



f 1 da a 1 (1 - a) r (a - a 2 ) n ^-f{ay + (1 - a)x) 
Jo da 



-(l + n) [l-l,r\n]dz f{z) + (r + n) [I, r - 1 1 n] dz f(z) . (2.68) 

J X J X 

In order to see the main difficulty, consider the example of two nested line integrals with 
two tangential derivatives 

(y - x) j {y-x) k f V [0, 1 | 0] dz t V™ (z t ) f [0, 1 | 0] dz 2 d jk V& (z 2 ) (2.69) 

J X J 21 

= (y - xy f [0,0 0] dZ! V^( Zl ) (y - Zl ) k At), 1 | 0] dz 2 d jk V^(z 2 ) 

J X J Z\ 

= -(y-x) j [ V dz! vWfa) djVWfa) (2.70) 

J x 

+(y - x) j [ V dzi V^{ Zl ) [ V dz 2 djV^iz^ . (2.71) 



Although the line integrals in ( 2.69| ) satisfy the conditions of Theorem |2.3| , the expression 
cannot be transformed into the required form ( |2.65 ). Namely in ( 2.70) ), we cannot eliminate 
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the remaining tangential derivative (because partial integration would yield a term {y — 
x)i djV^(zi)). In ( 2.71| ), on the other hand, we can successfully perform a second partial 
integration 

CT = f V [0,-l\0}dz 1 V^(zi)(V^(y)-V^(zi)) , 



but then the second parameter in the bracket [., . | .] becomes negative. More generally, 
we must ensure that the boundary terms contain no tangential derivatives, and that the 
parameters l a ,r a , and n a stay positive in the construction. 

Since the chirality of the potentials is not affected by the partial integrations, it is 
obvious that the rules (i) and (ii) of Theorem |2.3| will remain valid. For simplicity in 
notation, we will omit the indices c a in the following. 



First of all, we split up the factor (y — x) 1 in ( 2.32j ) in the form (y — x) = (y — x) K (y 



x) L , where L are those tensor indices which are contracted with the partial derivatives 
1, . . . , k. Setting 6 = 1 and zq = x, the first line integral in ( 2.32j ) can be written 



as 



d; a , a 



■■■ (y-z b ^) L I" [h,r b \n b ]dz b dilU%V^{z b ) ■■■ . (2.72) 

J Zb-l 

We rewrite the tangential derivatives in this line integral as derivatives in the integration 
variable, 



a 



d_ 

da 



d£>n?v£\z b ) 



with |L| 



= ••• (y- z^-x) / daa C 1 
Jo 

\N\ + q and I = l b + n b , r = r b + n b . Lemma gives the bounds 



I > 



1 



r > q+\N\ 



(2.73) 



(2.74) 



More generally, we use ( |2.73| ),( |2?74 ) as our induction hypothesis, whereby the left factor 
'■ • •' stands for all previous line integrals (which contain no tangential derivatives), and 
the right factor '■ • •' stands for subsequent line integrals. The tensor indices of the factor 
(y — z a -\) N must all be contracted with the partial derivatives d z a for a > b and thus 
give tangential derivatives in the subsequent line integrals. The induction step is to show 
that all the a-derivatives in ( |2.73[) can be eliminated, and that we can write the resulting 
expressions again in the form fl2,73j ),( t2.74 ) with b replaced by 6 + 1. Under the assumption 
that this induction step holds, we can eliminate all tangential derivatives in k steps. The 
resulting expressions are very similar to ( 2.65 ), ( 2. 66| ). The only difference is that the 



derivatives d Ka and n Pa in the resulting expressions are differential operators acting on all 
the following factors V^ a \ V^ a+l \. . . ; in ( 2.66j ), on the other hand, the partial derivatives 
act only on the adjacent potential V^ a \ In order to bring the resulting expressions into 
the required form, we finally carry out all the derivatives with the Leibniz rule and the 
chain rule (|2.49[) . 



For the proof of the induction step, we integrate in ( 2.73 ) q times by parts (if q is zero, 
we can skip the partial integrations; our expression is then of the form ( 2.76| )). Since the 
powers of the factors a and (1 — a) are decreased at most by one in each partial integration 
step, ( p.74| ) implies that the boundary values vanish unless in the last step for a = 0. We 
thus obtain a sum of terms of the form 



• • • (V ~ Zb-lf f [l,r\n = 0] dz b dg nil v 

J Zb-l 



and 



(2.75) 

with / > 0, r > \N\. (2.76) 
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In ( 2.76 ), we iteratively use the relation 



(y — x) 3 \ [I, r | n\ dz 



[l,r — 1 | 7i] (y — z) J 



(2.77) 



to bring all factors (y — to the right, 



(2.76) 



■f [l,r\n = 0]dz b (y- z b ) N d«» nil V J? 

Jz b _ 1 



with l,r > 0. (2.78) 



It might be convenient to reorganize the polynomials in this line integral with the identity 

rv 



[I, r | n] dz ■ 



rv 

/ [I - l,r - 1 [ n + 1] dz- 

J X 



but this is not relevant for the statement of the theorem. 

In both cases (2.75) and ( 2.78; ), we have an expression of the form 



(y - z b f (% □» V Jb 



(2.79) 



where the first factor '• • •' stands for line integrals without tangential derivatives, and 
where none of the factors (y — z b ) are contracted with d^ b b . Applying the "inverse Leibniz 
rules" 



(y 



d 
dx k 



d 

dx k 



{y - x) j + 5{ 



(y - x)j u x 



U x (y - x)j + 2^- 



(2.80) 
(2.81) 



we iteratively commute all factors (y — z b ) in ( 2.79 ) to the right. This gives a sum of 
expressions of the form 



(2.82) 



where the factors (y — z b ) are all contracted with the partial derivatives d* a , a = 6+1, ... ,k. 
The Leibniz rules may have annihilated some factors (y—z b ) (i.e., \L\ might be smaller than 
\N\); in this case, the parameters t a , a = b + 1, . . . , k have decreased. As a consequence, 
the inequalities of Lemma 2A are still valid for all expressions ( |2.82| ). If we write fl2.82p 
in the form ( 2.72|) with b replaced by b + 1, we can thus split up the tangential derivatives 
in the form ( 2.73| ),( p.74 ), which concludes the proof of the induction step. 

It remains to derive equation ( 2.6 7| ) : Note that each partial integration decreases both 
the number of factors (y — z a ^\) and the total number of partial derivatives by one. If we 
carry out the remaining derivatives with the Leibniz rule (in the last step of the proof), 
this does not change the total order J^o=i l-^al + 2p a of the derivatives. Therefore relation 
( |2.34| ) in Theorem |2.3| transforms into ( 2.67| ) . ■ 



2.2 Resummation, Reduction to the Phase-Free Contribution 

In the previous subsection, we gave a procedure for performing the light-cone expansion of 
every summand of the perturbation series ( 2,12j ) . In order to obtain formulas for the light- 
cone expansion of s, we must sum up the light-cone expansions of all Feynman diagrams. 



23 



Collecting the contributions of all Feynman diagrams gives, to every order on the light 
cone, an infinite number of terms. To get control over all these terms, we shall reorder the 
sums and partly carry them out. In the end, the light-cone expansion for s(x, y) will, to any 
order on the light cone, consist of only a finite number of summands. This rearrangement 
and simplification of the sums is called resummation of the light-cone expansion. 

In order to get a first impression of what needs to be done, we consider the leading 
singularity on the light cone (more precisely, we neglect all terms of the order 0((y — 
x)~ 2 )). This corresponds to taking only the contributions with h = — 1 in Theorem |2.3[ 
According to the bound fl2.36| ), the multi- indices I a and the parameters p a must all vanish. 
Furthermore, equation ( [2.34[ ) yields that 1 1\ = k + 1. The only possibility to satisfy the 
rules (a) and (b) is to contract one factor (y — x) with each potential V^, a = 1, . . . , k, 
and the remaining factor (y — x) with a 7-matrix. Thus the potentials must all be 
chiral potentials A L / R . According to the rules (a) and (b), the chirality of all potentials 
must coincide with the chirality of the projector Xc- We conclude that the leading order 
of ((—sB) k s)(x,y) on the light cone consists of a sum of expressions of the form 

ry 

XcC (y- x) h ■■■{y- x) jk I [h,n \ m] dz x A£{z\) 

J X 

ry 



[k,n I n k ] dz k A J c k (z k ) (y - x) j jj 5 ( x) (x, y) 



with c = L or c = R. Thus Theorem |2.3| makes a precise statement on the mathematical 
structure of all the contributions to the light-cone expansion. However, it does not give 
information about the values of the parameters C, l a , r a , and n a . In order to see more 
precisely how the leading order on the light cone looks like, we perform the light-cone 



expansion directly with Lemma 2.2, To first order in the external potential, we obtain 



-s B s)(x,y) 
£2 



(2.11), (2.25) 1 



- & / d 4 zS^(x,z) (B{z) (y - z) k lk ) S^(z,y) 



1 

2 j x 

1 ry 

llx 

1 ry 

2J X 



[0,-1 I 0] dzB(z) (y-z) k lk S^(x,y) 
dz (y - xYjj B(z) (y - x) k j k S^ix, y) 
dz (xl A{{z) + X r A j R {z)) (y - x)j (y - x) k lk S^ l \x,y) 
i(y-x)j f ' dz(x L A 3 L {z)+XRA R {z)) s(x,y) , (2.83) 



where 'x' denotes that we only take the leading order on the light cone. Since ( 2.83| ) is a 
product of a smooth function with a single factor s(x, y), this formula can be immediately 
iterated. We obtain for the left and right handed component of the k th order Feynman 
diagram 



Xc ((-« B) k s)(x,y) x xc (-if (V ~ x) h ■ ■ ■ (y - x) jk / [0, k - 1 | 0] dz x A£{ Zl ) 



[0, fe — 2 I 0] dz 2 Ai 2 (z 2 ) 
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[0, I 0] dz k Ai"(z k ) s(x,y) 



(2. 



Zfc-l 



(Notice that, according to Lemma |2.l| , the higher order terms which were neglected in 
( |2.83| ) also give contributions of the order 0((y — x)~ 2 ) in the iteration. Therefore it is 
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really sufficient to take only the leading contribution on the light cone in every step.) The 
line integrals in ( 2.84] ) are particularly simple. Namely, they are the k th order contributions 
to the familiar Dyson series. From this we conclude that, to leading order on the light 
cone, the sum over all Feynman diagrams converges absolutely. We can carry out the sum 
and obtain 



Xcs{x,y) x XcPexp^-i^ A 3 C \ ay+[1 „ a)x (y - x)j daj s(x,y) , 



(2.85) 



where Pexp is the usual ordered exponential. For completeness, we give its definition: 

Def. 2.6 For a smooth one-parameter family of matrices F{a), a € IR, the ordered expo- 
nential Pexp(J F(a) da) is given by the Dyson series 



Pexp 



(/>< 



a) da 



1 + f dt F(t ) dt + [ dt F(t ) f dt x F{t 

J a J a J to 

+ f b dt F(t ) [ dh F(ti) f b dt 2 F(t 

J a J to Jti 



b 
to 

2) + 



As is verified by a direct calculation, the ordered exponential is a solution of the ordinary 
differential equation 

-^-Pexp f J b F(a) daj = -F(a) Pexp f J* F(a) daj (2.86) 



with the boundary conditions 



Pexp (^J b F(a) da^j = 1 . 



(2.87) 



Because of the uniqueness of the solutions of ordinary differential equations, one can 
alternatively take ( |2. 86 ), ( 2.87 ) as the definition for the ordered exponential. 



For the ordered exponential in (2.85), we also use the shorter notations 



Pexp \-i j Al(y- x)j 



and 



Pe 



-i /" Al (y-x)j 



(21 



Notice that ( |2.88| ) is a unitary (/ x /)-matrix which depends only on the chiral potentials 
along the line segment xy. Its inverse isQ 



Pexp ( -i j A{{y- x)^] = Pexp ( -i J A 3 C (x - y)j 



(2.90) 



If the chiral potentials are commutative, i.e. [A c (x) , A c (y)] = 0, then the ordered expo- 
nential coincides with the ordinary exponential (this is e.g. the case if one considers the 
system / = 1 of only one Dirac sea). For the ordered exponential ( [2.88 ) along a line 



x We mention for clarity that this is in general not the same as 



Pexp [if Ai(y- x). 



(2.89) 



because the exponentials in (2.90) and (2.89) are ordered in opposite directions 
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segment in Minkowski space, the differential equation ( |2,86| ) can be written with partial 
derivatives as 

(y- x )k^L_p e -i£M(v-x) 3 = i( y _ x ) k A k c {x) Pe-*/"^^-^ . (2.91) 

We conclude that, to leading order on the light cone, the special form of the formulas 
of the light-cone expansion allows us to immediately carry out the sum over all Feynman 
diagrams. Unfortunately, the situation to higher order on the light cone is more difficult, 
because the combinatorics of the partial derivatives and of the tensor contractions becomes 
very complicated. We cannot expect that the sum over all Feynman diagrams can then 
still be written in a simple, closed form. Nevertheless, ordered exponentials over the 
chiral potentials should be helpful. More precisely, it is promising to write the light-cone 
expansion with line integrals which contain intermediate ordered exponentials, like e.g. 
the line integral 

V dz Pe-*£^ (p4 L (z)) Pe-'M {y - z)k (2.92) 

(expressions of this form are also suggested in view of the behavior of the fermionic pro- 
jector under local gauge transformations of the external potential). Our basic idea is to 
arrange the contributions to the light-cone expansion of Theorem to any given order on 
the light cone in such a way that all infinite sums (which arise from the fact that we have 
an infinite number of Feynman diagrams) can be carried out giving ordered exponentials. 
We want to end up with a finite number of terms which are of the form ( [2.65 ) with the 



only exception that the nested line integrals contain, similar to ( 2.92j ), additional ordered 
exponentials. 

Before we can make this idea mathematically precise, we must clarify the connection 
between the line integrals in ( |2.65| ) and the line integrals with intermediate ordered expo- 
nentials. For this, we consider the example ( p. 92 ). If we expand the ordered exponentials 



in a Dyson series and reparametrize the integrals, ( 2.92j ) goes over into an infinite sum of 



nested line integrals of the form as in ( p. 65 ), more precisely 



oo 



(2.92) = V [ V [h,r 1 \n 1 }dz 1 A J i(z 1 )(y-x) jl [ V [l 2 ,r 2 \n 2 }dz 2 A J i(z 2 )(y 



p,q=Q 
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• • • [ [l p , r v | n p ] dzp A j £(z p ) (y - x)j p j [l,r\ n] dz h4l( z ) 

y ■ py 

[l P +l,r p+1 | n p+1 ] dz p+1 A 3 [ +1 (z p+1 ) (y - x) jp+1 / A J £ +2 (z p+2 ) (y - x) jp+2 



z 



'Zp+1 

[lp+l , r P+q \ n p+1 ] dz p+q A J [ +q (z p+q ) (y - x) jp+q . (2.93) 

' z p+q — l 

The formulas of our desired light-cone expansion must be such that, after expanding the 
ordered exponentials in this way, we get precisely all the contributions to the light-cone 



expansion of Theorem |2.5| . We point out that we can view the expansion ( 2.93 ) as a power 
series in the functions A 3 L j R (y — x) j . The leading contribution to this power series is simply 
the line integral 

rv 

dzn4 L (z) ; 



it is obtained from (2.92) by taking out the ordered exponentials there. In view of this 



example, we can hope to get the contributions to our desired expansion without the ordered 
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exponentials (i.e. with all the ordered exponentials removed from the formulas) by picking 



those contributions to the light-cone expansion of Theorem 2.5 which contain no factors 



A 3 L , R (y — x)j. We take these contributions as the starting point for our construction. 
Def. 2.7 A contribution 63[ ), [2. 6(\ ) to the light-cone expansion of Theorem |g, <5| is called 



phase- free if all the tangential potentials Vj are differentiated, i.e. 



\K a \ + 2 Pa >0 



whenever 



J a is contracted with (y — x 



,K 



To leading order on the light cone, the only phase-free contribution is the free Green's 
function s (namely, according to (2.84), the contributions with k > 1 all contain factors 
■^l/r (y~ x )j)- The restriction to the phase- free contribution also simplifies the situation in 
the general case. Namely, the following proposition shows that the phase- free contributions 
of the higher order Feynman diagrams involve higher mass-derivatives of the Green's 
functions. 



Proposition 2.8 For every phase- free contribution (2.61) to the light-cone expansion of 
the k th order Feynman diagram ((—sB) k s)(x,y), the parameter h satisfies the bound 



h > -1 + 



k + 1 



(2.94) 



where [.] denotes the Gaufi bracket. 



Proof: Consider a phase-free contribution ( p.65|) of the k til order Feynman diagram. 
According to the rules for the possible contractions of the tensor indices, only one factor 
(y — x) may be contracted with <y J ; the remaining \K\ — 1 factors (y — x) must be contracted 
with the vj a \ . Thus at least \K\ — 1 potentials are tangential and must (according to 
Def. be differentiated. This gives the inequality 



\K\-l < J2(\ K a\ + 2Pa) ■ 
a=l 

We substitute this bound into ( |2. 67 ) and obtain 2h > — 2 + k; this is equivalent to ( |2.94| ). 



According to the explicit formula Q2.5p , the higher mass-derivatives of the Green's functions 
are of higher order on the light cone. More precisely, ( 2.94| ) and ( |2.6D yield that the phase- 
free contribution to the light-cone expansion of the k th order Feynman diagram is of the 
order 

0{(y-xf 9 ) with g = -2+ . (2.95) 



This means that, to every order on the light cone, only a finite number of Feynman 
diagrams contribute. As a consequence, there are, to every order on the light cone, only 
a finite number of phase-free terms. 

The phase-free contributions are useful because our desired light-cone expansion is 
obtained from them by inserting ordered exponentials into the line integrals. We do this 
"by hand," according to simple rules. 



Def. 2.9 For every phase-free contribution (2.61 ) to the light-cone expansion of Theo- 
rem we introduce a corresponding phase-inserted contribution. It is constructed 
according to the following rules: 
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(I) We insert one ordered exponential into each line integral and one ordered exponential 

at the very end. More precisely, the phase-inserted contribution has the form 

Xc C(y- x) K W {0) (x) ! V [h,Ti | m] dzi Pe^XT 1 A " ^~ x ^ w^{z{) 

Jx 

x f y [l 2 ,r 2 \n 2 ]dz 2 Pe- l ^ A ^ iZ2 ~ Z1 ^ W^\z 2 ) ••• 

J Z\ 

x [l a ,r a \n a ]dz a Pe W w (z a ) 

x p e - i /Idb- 2 »^ +1 ^ S^ h \x,y) . (2.96) 

(II) The chirality cp, (3 = 1, . . . , a + 1 o/ i/ie potentials in the ordered exponentials is 

determined by the number of dynamic mass matrices in the factors ; namely 

C/?_i and cr \ co ^ nc ^ e I jt w(fi—i) con t a i ns an ) even [ num \, er f factors Y, 
p p are opposite I ) odd I 



where cq := c zs t/ie chirality of the projector Xc in ( 2. 96 ). 



To illustrate these insertion rules, we consider the example of two nested line integrals 



XL 



d Zl (p4 L )(zi) / dz 2 mY L (z 2 ) SW(x,y) . 



The corresponding phase-inserted contribution is 

xl [ y d Zl Pe^r^-^Mn^X*!) 

^^ 2 Pe- l/ "^ (22 - 2l)fc m y L ( Z2 )P e - 1 ^^^ S W(x,y) . 



Zl 



Theorem 2.10 The light-cone expansion of the Green's function s(x,y) coincides with 
the sum of all phase-inserted contributions. 

Proof: A possible method for the proof would be to rearrange all the contributions to 



the light-cone expansion of Theorem |2.5| until recovering the Dyson series of the ordered 
exponentials in ( |2,96| ). However, this method has the disadvantage of being technically 
complicated. It is more elegant to use a particular form of "local gauge invariance" of 
the Green's function for the proof: First we will, for given x and y, locally transform 
the spinors. The transformation will be such that the light-cone expansion for the trans- 
formed Green's function s(x, y) consists precisely of all phase-free contributions. Using the 
transformation law of the Green's function, we then show that the light-cone expansion of 
s(x, y) is obtained from that of s(x, y) by inserting unitary matrices into the line integrals. 
Finally, we prove that these unitary matrices coincide with the ordered exponentials in 
Def. fO]. 

In preparation, we consider the transformation law of the Dirac operator and the 
Green's function under generalized local phase transformations of the spinors: We let Ul 
and Ur be two unitary matrix fields acting on the Dirac sea index, 

U L = (U L (x) l m \ m=1 _ f , U R = {U R {x) l m \ m=l ,...j with U L U* L = 1 = U R U* r , 
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and transform the wave functions Vl/ = {^f{x))f = i , '" , j according to 

*(a?) -> = C/"(a;)*(s) with Cf(x) = xlU l (x) + xrUr(x) • 

Thus t/x, and ?7r independently transform the left and right handed component of the 
wave functions, respectively. Notice that the transformation U is not unitary with respect 
to our scalar product (|1.4j) , because 

V := U- 1 = xl Uj} + X R Ur 1 but 
U* = 7° 7° = XR Ul 1 + XL Ur 1 . 

Therefore we must in the following carefully distinguish between U, U* and their inverses 
V, V*. As an immediate consequence of the Dirac equation (i@ + B — m)^ = 0, the 
transformed wave functions \l/ satisfy the equation 

V*(i0 + B)V% = . 

A short computation yields for the transformed Dirac operator 

V*($+B)V = i$+B 

with 

B = Xl(4r-™Yr) + X R{h-mY L ) , (2.97) 
whereby Ajm and Y L / R are the potentials 

*U = U L/R A J L/R UllR + iU L/ R{dWZ?R) (2.98) 
Yl/r = U L/R Y L/R U^ L . (2.99) 

We denote the advanced and retarded Green's functions of the transformed Dirac operator 
+ B by s. They satisfy the equation 

(ift. + B(x))s(x,y) = 5 A {x-y) . (2.100) 

Since we can view B as the perturbation of the Dirac operator, the Green's function s is, 
in analogy to ( 2.12j ), given by the perturbation series 

oo 

s = ^(-sB)"s . (2.101) 

n=0 

The important point for the following is that the Green's functions s and s are related to 
each other by the local transformation 

s(x,y) = U(x)s(x,y)U(y)* . (2.102) 



This is verified as follows: The right side of (2.102) also satisfies the defining equation 
( |2.100| ) of the Green's functions; namely 

{i$ x + B(x))U(x)s{x,y)U(y)* = V(x)* {i$ x + B(x)) V(x) U(x) §{x, y) U(y)* 
= V{x)*{i$ x + B{x))~s(x,y)TJ{y)* = V {x)* 8\x - y) U (y)* 
= V(x)* U{x)* 5 4 {x - y) = 5 A {x-y) . 
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Furthermore, the supports of both sides of ( |2.102| ) lie (depending on whether we consider 
the advanced or retarded Green's functions) either in the upper or in the lower light cone. 
A uniqueness argument for the solutions of hyperbolic differential equations yields that 
both sides of ( 2.102 ) must coincide. 

Next, we specify the unitary transformations Ul and Ur: We fix the points x and y. 
For any point z on the line segment xy, we chose Ul/r(z) as 



U l /r(z) = Pexp \-i I A° L/R (z - x) 
The differential equation ( 2.91| ) yields that 

(y-xyU c (z) (djUciz)- 1 ) 



(2.103) 



i /" A* (x-z) k 



Pe^I^c (*-*)* i(jJ _ x ^ A ^ Pe -iJ"Ai (x-z) k 
i(y - x)j U c (z) Al(z) U^z)- 1 . 



Substituting into ( 2.98 ) gives 

A J L/R (z) (y - x)j 



for z € xy. 



(2.104) 



Thus our choice of Ul and Ur makes the potentials Al(z) and Ar{z), z E xy, orthogonal 
to (y — x). Notice, however, that since Y l /r in ( 2.99|) is arbitrary and independent of Ul/r, 
( |2.103| ) gives no constraints for the dynamic mass matrices Yl/r- 

We point out that we did not specify Ul/r{z) outside the line segment z € xy; the 
unitary transformation Ul/r may be arbitrary there. This means that also A L / R is un- 
determined outside the line segment xy. Especially, all the non-tangential derivatives of 
A l /r(z), z G xy, are undetermined. However, equation (2.103) does give constraints for 
the tangential derivatives. For example, differentiating (2.104) in the direction (y — x) 
yields 

(y - x) j {y - x) k djA k L/R (z) = for z € xy . (2.105) 



We now consider the perturbation expansion ( 2.101 ). The light-cone expansion of all 
Feynman diagrams according to Theorem 2.5 gives a sum of terms of the form 



XcC (y- x) 



K 



W {0) (x) 

V 



[h,r 1 \n 1 ]dz 1 wWfa) 



[l a ,r a | n a ] dz a W^(z a ) 7 J S^ h \x,y) 
where the factors W(@) are given by 



[h,r 2 \n 2 ] dz 2 W ( - 2) (z 2 ) 

(2.106) 



WW = (d Ka f i U Pa f i V {af3) 



(2.107) 



Because of condition ( |2.104| ), all the contributions which are not phase- free vanish. Fur- 
thermore, according to Theorem 2.5 , the contributions ( 2.106 ), ( 2.107| ) contain no tan- 
gential derivatives. Clearly, the derivatives in these formulas may have a component in 
direction of {y — x); but the contribution of the derivatives transversal to (y — x) uniquely 
determines the form of each derivative term. Therefore all the phase-free contributions of 
the form ( 2.106|) , (|2.107|) are independent in the sense that we have no algebraic relations 
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between them. We conclude that, as long as the potentials A l /r and Y l /r are only specified 
by ( pT98| ), ( p9|) , and ( plQ3|) , the light-cone expansion pl06| ), pl07| ) consists precisely 
of all phase- free contributions. 

Next, we exploit the local transformation law ( |2.102| ) of the Green's functions: We 
solve this equation for s, 



s(x,y) = V(x) s(x,y)V(yY 



(2.108) 



The transformation Ul/r does not enter on the left side of this equation. Thus the right 
side of ( |2.108; ) is also independent of Ul/r- Especially, we conclude that the light-cone 
expansion of s(x, y) must be independent of the derivatives of Ul/r along the line seg- 
ment xy. At first sight, this might seem inconsistent because the individual contributions 
( |2.106| ),( [2.107| ) do depend on the derivatives of Ul/r (this is obvious if one substitutes 
fl2.98p and ( |2.99| ) into ( |2.107 ) and carries out the derivatives with the Leibniz rule). The 
right way to understand the independence of s(x, y) on the derivatives of Ul/r is that all 
derivative terms of Ul/r cancel each other to every order on the light cone if the (finite) 
sum over all contributions ( 2. 106| ) to the light-cone expansion of s(x,y) is carried out. 
Since we will form the sum over all contributions to the light-cone expansion in the end, 
it suffices to consider only those contributions to the light-cone expansion which contain 
no derivatives of Ul/r- This means that we can substitute ( gjgp and ( gjgg) into ( p.l07|) , 
forget about the derivative term iU L /R(&>U£, R ) in ( |2.98| ), and pull the unitary transfor- 
mations Ul/r, U£, r out of the derivatives. In other words, we can replace (j2lQ7p , 
by 



W(P) = U Ca (d Ka m Pa PV 



d a . 



K, 



(2.109) 



with chiral indices c a ,d a = L/R. The light-cone expansion for s(x,y) consists precisely of 
the sum of all phase- free contributions of the form ( |2.106[) , ( |2.109| ) . 

The chiralities c a , d a of the unitary transformations Ul/r, Uy R in ( |2. 109| ) are deter- 
mined by the rules (i) and (ii) (in Theorem |2.3| ) and by (2.98) and ( p. 99 ). According to 
the rule (ii), the indices c a and c a +i coincide iff is a chiral potential. According to 
( |2.98|) and (|2.99| ), on the other hand, the indices c a and d a coincide iff V^ a > = A l /r- We 
conclude that the indices d a and c a+ \ always coincide. Thus all the intermediate factors 
Ud a U Ca+1 give the identity, and ( |2.109| ) simplifies to 



WW = Uc0 W W u-i . 



(2.110) 



Furthermore, the chiralities ca and da coincide if and only if contains an even number 
of dynamic mass matrices. 

Finally, we substitute the light-cone expansion ( 2.106| ), (2.110) for s(x,y) into ( |2,108| ). 
This gives for the light-cone expansion of s(x, y) a sum of expressions of the form 



Xc C (y - x) K U-\x) (U C0 W^U^)(x) / [l 1 ,r 1 \n 1 ]dz l (U cl W^U^ l )( Zl ) 

J X 

[L,r a I n a ] dz a (U Ca W^U^)(z a ) U Ca+1 (y) 7 J S^{x,y) , (2.111) 



v 

Z<x-1 



where the sum runs over all phase-free contributions of this type. Similar to the consid- 
erations before ( 2.110D , one sees that adjacent unitary transformations always have the 
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same chirality. Thus ( |2.111| ) 

can be simplified to 



Xc C(y- x) K W<®(x) [h,n\ m] dz x u ci {x)- 1 U Cl ( Zl ) 

J X 

[l a ,r a | n a ] dz a U^iza-iy 1 U Ca (z a ) W {0 \z a ) U Ca+1 (z a y l 
xU Ca+1 (y)l J sW(x,y) , 



y 



whereby the indices c a satisfy the rule (II) of Def. 2^9. According to ( |2.103| ), the factors 
U ( T 1 (.) U c {.) coincide with the ordered exponentials in ( 2.96 ), which concludes the proof. 



For clarity, we point out that all the constructions following Def. 2.9 are based on the light 



cone expansion of Theorem 2.5. It is essential for the statement of Theorem 2.10 that the 
phase-free contributions contain no tangential derivatives. If we had worked with the 



light-cone expansion of Theorem |2.3| (instead of Theorem 2.5), the light-cone expansion 
of s(x,y) would not have consisted of all the phase-free contributions to the light-cone 
expansion. For example, a line integral containing the tangential derivative ( |2. 105 ) would 
vanish, although it is phase-free. As a consequence of this problem, the whole construction 
would break down. 

The introduction of the phase-free and phase-inserted contributions has simplified the 
light-cone expansion of the Green's functions considerably: Assume that we want to per- 
form the light-cone expansion to some given order on the light cone. Then we first calculate 



the phase-free contribution to the light-cone expansion; according to Proposition 2.5, this 



gives only a finite number of terms. Using the rules of Def. 2.9, we can easily construct the 



corresponding phase-inserted contributions. According to Theorem |2.10| , this finite num- 
ber of phase-inserted contributions gives precisely the light-cone expansion of the Green's 
functions to the desired order on the light cone. This procedure is called the reduction to 
the phase-free contribution. 



2.3 Calculation of the Phase- Free Contribution 

According to the reduction to the phase-free contribution, it remains to calculate the 
phase-free contribution to any given order on the light cone. Although this is still a very 
complicated problem, we know from Proposition 2.8 that we only get a finite number of 



terms. This makes it possible to use a computer algebra program for the calculation. The 
author has developed the C++ program "class_commute" specifically for this problem. It 
generates explicit formulas for the phase-free contribution to any order on the light cone. 
We now outline how this program works, without entering implementation details^. 

All the objects occurring in the calculation (like integrals, partial derivatives, Lapla- 
cians, potentials, Dirac matrices, etc.) are described by different data structures (classes 
in C++). Formulas are built up as sums of lists of these data structures. The calculation 
is performed by manipulating the lists. More precisely, this works as follows: Each data 
structure carries an ordering number. At the beginning of the computation, the lists are 
disordered in the sense that their elements do not occur with increasing ordering numbers. 
For ordering the lists, the program iteratively commutes adjacent elements of a list. Each 

2 The commented source code of the program "class_commute" is available from the author on request. 
It is an extension of the program used in g| for the light-cone expansion to first order in the external 
potential. 
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commutation is performed by a function of the program which is specific to the particular 
pair of data structures; i.e., there is a function for commuting a partial derivative with a 
potential, a function for commuting two Dirac matrices, etc. (this is easily implemented 
in C++ using virtual class functions). The data structures and commutation rules are 
designed in such a way that, after the ordering process has come to an end, the lists consist 
of the desired formulas of the light-cone expansion. 

The main advantage of this implementation with commutation rules is that the pro- 
grammer must only think of the calculation on a "local" scale by telling the computer the 
rules for commuting a given pair of data structures. Furthermore, this gives a convenient 
segmentation of the computer program into small, independent parts, which can be writ- 
ten and debugged separately. As soon as all commutation rules are specified correctly, the 
program can perform the whole calculation "globally" by recursively applying the com- 
mutation rules. Compared to standard computer algebra packages like e.g. Mathematica 
or Maple, this concept is very flexible and efficient. 

For the more complex manipulations like partial integrations and the handling of tensor 
indices, the program uses a so-called "message pipe," through which the objects in the 
formulas can pass information to other objects (this is again implemented with virtual class 
functions). In this way, the elements of the lists can exchange data and give commands 
to each other. This allows a convenient coordination of the formula manipulations. 

The calculation rules of the program "class_commute" are very similar to the construc- 
tion steps described in the proofs of Theorem 2.3 and Theorem 2.5. The only difference 
is that, according to our implementation as commutation rules, the calculation does not 
follow the same strict and clear order as in Subsection 2A_. Basically, one may think of the 
construction steps of Theorem 2.5 and Theorem 2.5 as being performed simultaneously in 
a disordered way, whereby the program ensures that all rules are applied consistently. 

Some formulas generated by the program "class_commute" are compiled in the ap- 
pendix; they give a picture of the leading singularities of s(x, y) on the light cone. We 
remark that "class_commute" was also a valuable tool for finding and checking the com- 



binatorial results of Theorem 12.31, Lemma 2.4, and Theorem 2.5 



3 The Light- Cone Expansion of the Dirac Sea 

In this section, we shall perform the light-cone expansion of the fermionic projector as 
defined in [ffl. In Subsection |3.1| , we establish a formal analogy between the light-cone 
expansions of the Dirac sea and of the Green's functions. This residual argument allows 
us to use the results of the previous section also for the fermionic projector. However, 
the analogy between the Dirac sea and the Green's functions cannot be extended beyond 
a purely formal level. The basic reason is that, in contrast to the Green's functions, the 



Dirac sea is a non-causal object. This is developed in detail in Subsections [O] and 3.3. 

We point out that, in this section, we do not work with the dynamic mass matrices 
Yjj r (x), ( |2.8D . The reason is that, for the Dirac sea, the regularity conditions of Lemma 
|0] are necessary for the contributions to the perturbation expansion to be well-defined. 
Working with the dynamic mass matrices, however, implies that we consider the potential 
B, ( |2.10| ), as the perturbation of the Dirac operator; but B does in general not go to zero 
at infinity. For our notation, the reader is referred to fl|. 
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3.1 The Residual Argument 

We begin by describing how the light-cone expansion of the Green's functions can be 
understood in momentum space. Apart from giving a different point of view, this allows us 
to get a connection to the light-cone expansion of the Dirac sea. For clarity, we begin with 
the special case mY = of zero fermion mass. This case is particularly simple because then 
B = £>, so that the perturbation expansions Q2.1| ) and ( 2.12| ) coincide. This is sufficient to 
explain the basic construction; the extension to mY ^ will later be accomplished by a 
general argument. Furthermore, we only consider the advanced Green's function; for the 
retarded Green's function, the calculation is analogous. 

Suppose that we want to perform the light-cone expansion of the k th order contribution 
to the perturbation series (2.1)=( p.l2 ). We write the contribution as a multiple Fourier 
integral, 

((-S V Bf S V)( X ,y) 

d* P f d\ t f dUk Afi V (;p;(7i) __ ^ qk)e - i{p+qi + ... + qk)x + lp y ^ (31) 



(2tt) 4 7 (2tt) 4 J (2vr) 4 
where the distribution As v (p; q%, . . . , q k ) is the Feynman diagram in momentum space, 

As v (p; qi,...,q k ) = (-l) k s v (p + q x + • • • + q k ) B(q k ) s v (p + gi + ■ • • + q k - X ) B(q k -i) 

••• B{q 2 )s v ( P + q 1 )B(q 1 )s v {p) (3.2) 

(-B denotes the Fourier transform of the potential B, and s v (p) is the multiplication 
operator in momentum space). For the arguments of the Green's functions, we introduce 
the abbreviation 

p := p and pi :=p + qx-\ \- qi , 1 < I < k. (3.3) 

Substituting the explicit formulas ( ^.3[) and ( 2 . 1 5|) into (^^ ), we obtain 
As v (p; qi ,...,q k ) = (-l) k i k B{ Slk )i k _ x ■ ■ ■ i x B( Ql ) p* 

x lim 1 1 1 (3 4) 

o< £ ^o ( Pk ) 2 - iep° k (p k -i) 2 - iepl_i (p ) 2 - iep% 

We now expand the Klein-Gordon Green's functions in ( |3.4| ) with respect to the momenta 
Pi — p. If we expand the terms iep® with a geometric series, 

1 ~ (ie(p°-p°)) n 



E 



{pi) 2 - iep\ ^ Q {{pif - iep°) 1+n 

all contributions with n > 1 contain factors e and vanish in the limit e — > 0. Therefore 
we must only expand with respect to the parameters ({pi) 2 —p 2 )- This gives, again with 
geometric series, 

A S v (p; qi ,...,q k ) = (-l) k ^ k B(q k )^ k _ 1 ■ ■ ■ i x B{q x ) p* Q 



x 



E (p 2 -plr k ■■■ (p 2 -pi) ni J™ 



0<e^0 (p 2 - i £ pOY+k+ ni +-+n k 
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Rewriting the negative power of (p 2 — iep°) as a mass-derivative, 

1 



(jj2 — ^gpOU+fe+niH hn fc 

1 / f£ \ fc+niH — l_rafc 1 



{k + n\ + ■ ■ ■ + rik)\ \da) p 2 — a — iep° i o _ 



(3.5) 



we obtain a formula containing only one Green's function. Namely, using the notation 

(H), 

As v (p; qi,... ,q k ) = (-if j> k B{q k ) p 1 ^ ■■■j 1 B(q 1 )j Q 

OO -I 

x £ on z — zrnr (p 2 -^) nfc • • • (p 2 ~p 2 i) ni s v(fc+ni+ - +nfc) (p) . (3.6) 

m,...,n fc =o ^ + 1 h nfc -'- 

This is the basic equation for the light-cone expansion of the Green's functions in mo- 
mentum space. Similar to the light-cone expansion of the previous section, (|3.6| ) involves 
the differentiated Green's functions S v ^ . It remains to transform the polynomials in the 
momenta po, ■ ■ ■ ,Pk until getting a connection to the nested line integrals of, say, Theorem 
Substituting (3.3), we rewrite (3J5) in terms of the momenta p, q±, . . . , q k and multiply 
out. Furthermore, we simplify the Dirac matrices with the anti-commutation rules ( 2.530 . 
This gives for (3.6) a sum of terms of the form 

XeCYgf Vj"lM)--^.hlM)P L SV(h) (p) (h>[\L\/2}), (3.7) 

where the tensor indices of the multi-indices /, Ii, J;, and L are contracted with each 
other (similar to the notation of Theorem the factors Vj c stand for the individual 
potentials of B). If tensor indices of the power p L are contracted with each other, we can 
iteratively eliminate the corresponding factors p 2 with the rule ( p. 14 ), more precisely 

p 2 S v(h \p) = frS v(ft-1) (p) (h>l). (3.8) 

Thus we can arrange that the tensor indices of p L in ( |3.7D are all contracted with tensor 
indices of the factors 7 7 , q\\ or VjX.- By iteratively applying the differentiation rule 
( |2.16| ), we can now rewrite the power p L in ( |3.7D with p-derivatives, e.g. 

P,PkS-m(p) = -I Sj ^ t svw w = -iJ t o> J svWo,)) + i Sjt svw w 

In this way, we obtain for As v (p; q\, . . . , q k ) a sum of terms of the form 

Xc C 7 J ^- flf 1 V^ Ck (q k ) • • • V^ cl {q x ) fif S v W(p) , (3.9) 

whereby no tensor indices of the derivatives dp are contracted with each other. We 
substitute these terms into ( |Q| ) and transform them to position space. Integrating the 
derivatives dp by parts gives factors (y — x) K . The factors qj l , on the other hand, can be 
written as partial derivatives d Jl acting on the potentials V®. More precisely, the term 
fl3.9|) gives after substitution into (|3.1|) the contribution 

XcCi \h\+-+\h\ { _i)\K\ y (d I ^l k (x))---(d I wil ) iCi (x)) (y-x) K S v( - h \x,y) , (3.10) 
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where the tensor indices of the factor (y — x) are all contracted with tensor indices of 
the multi-indices /, Ii, or J/. The Feynman diagram ((—sB) k s)(x,y) coincides with the 
sum of all these contributions. 

This expansion has much similarity with the light-cone expansion of Theorem [2.3| . 
Namely, if one expands the nested line integrals in (12. 321) in a Taylor series around x, 



one gets precisely the expansion into terms of the form ( |3 . 1 0|) . Clearly, the light-cone 
expansion of Theorem ^ goes far beyond the expansion ( 3.1C| ), because the dependence on 



the external potential is described by non-local line integrals. Nevertheless, the expansion 
in momentum space (3.6) and subsequent Fourier transformation give an easy way of 



understanding in principle how the formulas of the light-cone expansion come about. We 
remark that, after going through the details of the combinatorics and rearranging the 
contributions ( 3.10| ), one can recover the Taylor series of the line integrals in ( |2.32| ). This 



gives an alternative method for proving Theorem 2.3. However, it is obvious that this 
becomes complicated and does not yield the most elegant approach (the reader interested 
in the details of this method is referred to ||], where a very similar technique is used for 
the light-cone expansion to first order in the external potential). 

Our next aim is to generalize the previous construction. Since we must, similar to 
( |3.5D , rewrite a product of Green's functions as the mass-derivative of a single Green's 
function, we can only expect the construction to work if all Green's functions in the 
product ( |3.2| ) are of the same type (e.g. the construction breaks down for a "mixed" 
operator product containing both advanced and retarded Green's functions). But we need 
not necessarily work with the advanced or retarded Green's functions. Instead, we can 
use Green's functions with a different position of the poles in the complex p°-plane: We 
consider the Green's functions 

s ± (p) = dS^ Jp) with S ± (p) = lim -= — (3.11) 

and again use the notation ( 2.13j ), 



o±(0 - ( jL\ <?± 

The perturbation expansion for these Dirac Green's functions is, similar to ( p.l| ) or ( |2,12| ), 
given by the formal series 

oo oo 

s+ := Y.(~ s+B ) ns+ > *~ : = J2(~ s ~ B T S ~ ■ ( 3 - 12 ) 

n=0 n=0 

The light-cone expansion in momentum space is performed exactly as for the advanced 
and retarded Green's functions; we obtain in analogy to ([O]) and ( |3.6|) the formula 

((- s ± B) k s±)(x,y) 

d A p f d 4 qi f d 4 q k ± 



(2vr) 4 J (2vr) 4 J (2tt) 



with 

,±r^. „ „ \ — I -t\k 



As ± {p;qi,...,q k ) = (-l) fc i> k B{q k ) i> k _ x ■■■$ 1 B{q 1 )$ Q 



oo 

X 

n\,...,n k 



OO -i 

y ( P 2 - p\) nk ■■■ (p 2 - P i) ni s ± 



36 



bmce S ± are Green's functions of the Klein-Gordon equation, they clearly also satisfy the 
identity (3.8). Furthermore, the differentiation rule ( p. 16 ) is also valid for S , namely 



9 ^\ P ) = hm 9 ( ' 



dpi \da ) o<e-+o dp> \p 2 — a z fie 



a=0 



Pi' lim t 2 ~ 2?3 ■ ^2 = -2P 1 S ± « +I \p) 
da J 0<s->o (p 2 - a =F is) 2 i o=0 



Therefore we can, exactly as in ( |3.9|) , rewrite the power p L with p-derivatives. Thus the 
expansion ( |3.10| ) is valid in the same way for the Green's functions s if one only replaces 
the index " v " in (glc|) by a± ". As explained before, the expansion (|3.1C| ) is obtained from 



the light-cone expansion of Theorem 2.3 by expanding the potentials around the space- 



time point x. Since the formulas of the light-cone expansion are uniquely determined by 



this Taylor expansion, we immediately conclude that the statement of Theorem 2.3 is also 
valid for the k th order contribution to the perturbation expansion (CT) if the factor m 
(IPD stands more generally for S + W or S ^ , respectively. This simple analogy between 
the formulas of the light-cone expansion of the Feynman diagrams (( — s v / A B) k s v / A ) and 
((— B) k s^), which is obtained by changing the position of the poles of the free Green's 
functions in momentum space, is called the residual argument. 

After these preparations, we come to the fermionic projector in the general case mY ^ 
0. We want to extend the light-cone expansion to an object P res being a perturbation 
of the free fermionic projector. Our method is to define P ies in such a way that it can 
be easily expressed in terms of the Green's functions s v , s A , s + , and s~ . The light-cone 
expansion of the Green's functions then immediately carries over to P res . We denote the 
lower mass shell by T a , i.e. in momentum space 

T a {q) = @(-q°) 5(q 2 - a) , (3.13) 

and set 

TV) - (s)' r -i~ ■ < 3 - 14 > 

Furthermore, we introduce, exactly as in pi, the series of operator products 



EH^ , b = £(-£s) fc £ , b> = 

fc=0 k=0 k=0 



Def. 3.1 The residual fermionic projector P res (x,y) is defined by 

P res (x,y) = \x(p™-~k){x,y) , (3.15) 
where the operators p res and k are given by the perturbation series 

oo 

-res = Y^{-iTt) 2 P b < p{bpfH > (3.16) 
(3=0 
oo 

k = j2(- i7T ) 2Pb<k ( bk ) 2Pb> ■ ( 3 - 17 ) 
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Proposition 3.2 (formal light-cone expansion of the residual fermionic projec- 
tor) The results of Section also apply to the residual fermionic projector. More precisely, 
the light-cone expansion of Theorem <L; holds for P res (x,y) if we replace by and 
multiply the formulas of the light-cone expansion from the left with the chiral asymmetry 
matrix X . According to Theorem 2.1, all tangential derivatives can be integrated by parts. 
With Def. 13 Def. |J, and Theorem \2Adi , the light -cone expansion can be reduced to the 
phase-free contribution. According to Proposition \2.S{ , the phase-free contribution consists, 
to every order ~ T^ 1 ) ( of only a finite number of terms. 

Proof: First of all, we must generalize the residual argument to the case mY ^ of 
massive fermions. According to ( |2.1| ) and ( 2.1 2| ) , there are two equivalent perturbation 
series for s v , 



J2(s v Bfs" 1 

k=0 

oo 

k=0 



(3.18) 
(3.19) 



In both perturbation series, each summand is a well-defined tempered distribution (this 
follows from the smoothness of B, B and from the causality of the perturbation expansion). 
In Section ^, we developed the light-cone expansion from the series in ( 3.19j ). But by 
arranging the contributions to this light-cone expansion in powers of the potential B, 
one also obtains formulas for the light-cone expansion of every Feynman diagram of the 
perturbation series ( 3.18; ). For the Green's functions s^, we have similar perturbation 
expansions 



k=0 

oo 

= s m=0 B) S m =0 

k=0 



(3.20) 
(3.21) 



Since the support of the distributions s ± (x,y) does not vanish outside the light cone, we 



now need the conditions of Lemma 1.1 on the decay of the potentials at infinity. According 
to our assumptions on B, each summand of the perturbation expansion ( p. 20 ) is a well- 
defined distribution. The potential B, however, does in general not decay at infinity; 
thus the Feynman diagrams of the perturbation expansion ( |3.21| ) are ill-defined. This 
is a problem, especially because in our above consideration, the residual argument was 
derived for the Feynman diagrams of the expansions ( j3l9D and ( ^2l| ). The solution to 
this problem is an approximation argument using the "causality" of the formulas of the 
light-cone expansion: We consider a smooth function ?7_r(x) which is equal to one inside 
the ball of radius R around the origin and vanishes outside the ball of radius 2R (in IR 4 
equipped with the standard Euclidean metric). Then the potential n^B has compact 



support and, according to Lemma 1.1, the Feynman diagrams 



k „± 



(3.22) 



are well-defined. We can apply the above residual argument for mY = 0; this yields formu- 
las of the light-cone expansion in terms of the potential {rjnB) and its partial derivatives. 
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Since the potential enters into the formulas of the light-cone expansion only along the 
convex line xy , we can, by taking the limit R — > oo, remove the cutoff function rjn from 
these formulas. This limiting process shows that the summands of the perturbation series 
in ( 3.2 1| ) make mathematical sense in terms of the light-cone expansion. By reordering 
the contributions, we immediately also get formulas for the light-cone expansion of the 
Feynman diagrams of the perturbation series ( [3.20 ), The analogy between the light-cone 
expansions of the Feynman diagrams of the perturbation series ( [3.18 ) and (3.20) finally 
yields the extension of the residual argument to a general mass matrix mY. 

Evaluating the poles in ( 2.15| ) and ( p . 1 2j ) in the complex p°-plane gives (using the 
formula limo< e ^o(( a; — ie) _1 — (x + is)^ 1 ) = 2mS(x)) the relations 



s + iir k , 
s + iir p , 



s 
s 



iir k 
iir p 



(3.23) 
(3.24) 



where s denotes as in 
functions, 



the arithmetic mean of the advanced and retarded Green's 



v + s A ) 



We substitute ( 3.23 ) and ( 3.24 ) into the perturbation series ( |2.lD , ( 3.20 ) and multiply out. 
After rearranging the sums, one sees that the series ( 3.16[ ) and ( |3.17| ) can be written as 



pves 



— (s+ - s~) 
2m y 



and 



1 
2~ki 



— (s v - s A ) 



(3.25) 



respectively (see JD, proof of Theorem 3.2] for the details of the combinatorics). According 
to the residual argument, all Green's functions have a light-cone expansion according to 
Theorem 2.3, By substituting into ( 3.25| ), this light-cone expansion immediately general- 
izes to p res and k. Using ( 3.15| ), we conclude that Theorem is also valid for P res after 



the replacement — * and multiplication with the chiral asymmetry matrix. Since 
the results of Theorem |2.5| , Proposition 2.8, and Theorem p.10 are obtained merely by 



manipulating and rearranging the formulas of the light-cone expansion, they also hold for 
the residual fermionic projector. ■ 



We point out that the argumentation in this subsection was only formal in the sense that 
we did not care about the convergence of the infinite sums. Also, the approximation argu- 
ment in the proof of Proposition 3^ requires a mathematical justification. Furthermore, 
the derivative ( |3,14 ) is ill-defined because, for a = 0, the mass shell degenerates to the 
cone {q 2 = 0, q° < 0}, which is not differentiable at q = 0. We postpone the mathematical 



analysis of these problems to Subsection 3.5. 



3.2 The Non-Causal High Energy Contribution 

Unfortunately, the residual fermionic projector P res of the previous subsection does not 
coincide with the fermionic projector P of Q, 

1 



P(x,y) 



X(p-k)(x,y) . 



(3.26) 



The difference is that, instead of the operator p Tes in the residual fermionic projector 
( |3.15| ), the fermionic projector ( |3.26 ) involves the operator p, which is formally given by 



P 



formally 



(3.27) 
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Using an operator calculus method, this formal definition is made mathematically precise 
in [0 in terms of a perturbation series for p. Similar to ( |3.16 ), this perturbation expansion 



consists of a sum of operator products. But the operator products are more complicated; 
they also contain operators k with some combinatorial factors (see [|]] for details). 

Before entering the mathematical analysis of the operator products, we point out 
that it is not just a matter of taste to take ( 3.26| ), and not ( 3,15| ), as the definition of the 



fermionic projector; only the definition ( 3.26| ) makes physical sense. This comes as follows: 



As explained in the operator k generalizes the splitting of the solutions of the Dirac 
equation into solutions of positive and negative frequency to the case with interaction. 
The "generalized positive and negative frequency solutions" are given by the eigenstates 
of k with positive and negative eigenvalue, respectively. The construction ( 3.27 ) , ( pT26| ) 



projects out all positive eigenstates of k; the operator |(p — k) consists precisely of all 



eigenstates of k with negative eigenvalue. The residual fermionic projector ( |3.15 ), however, 



consists of a mixture of positive and negative eigenstates of k, which is not a reasonable 
physical concept. 

We begin by giving the difference between the fermionic projector and the residual 
fermionic projector a name. 

Def. 3.3 The non-causal high energy contribution P he (x,y) to the fermionic pro- 
jector is given by 

P he (x,y) = P(x,y) - P res (x,y) . 

Clearly, this definition is only helpful if P he has some nice properties. The reason why 
the definition makes sense is that every contribution to the perturbation expansion of 
P^ c (x,y) is a smooth function. Thus the singular behavior of the fermionic projector on 
the light cone is completely described by the residual fermionic projector and its light-cone 
expansion, Proposition |3.2|. 



Theorem 3.4 The non-causal high energy contribution P he (x, y) is, to every order in 
perturbation theory, a smooth function in x and y. 

Proof: The perturbation series [Q, Theorem 4.1] defines p as a sum of operator products 
of the form 

CnBCn-xB ••■ BC , (3.28) 

where the factors C/ coincide with either k, p, or s. The number of factors k in these 
operator products is always even. If one replaces all factors k by p, one gets precisely the 
perturbation series for p Tes , ( 3.16| ) (this is verified using the details of the combinatorics in 



[yj ). Therefore we can convert the perturbation series for p into that for p res by iteratively 
replacing pairs of factors k in the operator products by factors p. Thus the difference 
(p— p res ) can, to every order in perturbation theory, be written as a finite sum of expressions 
of the form 

C n B---C b+1 B( P BC b ^---C a+1 Bp 

-kBC b ^---C a+1 Bk)BC a ^---BC , (3.29) 

where the factors C\ again stand for k, p, or s. Since P he = I X(p - p res ), it suffices to 
show that Q3.29 ) is a smooth function in position space. 



We first simplify our problem: Once we have shown that the bracket in (|3.29|) is smooth 



and bounded in position space, the additional multiplications to the very left and right 
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can be carried out by iteratively multiplying with B and forming the convolution with Ci, 
which again gives a smooth and bounded function in each step (notice that, according to 
the assumptions of Lemma LI , B decays sufficiently fast at infinity). Thus we must only 
consider the bracket in (3.29). We rewrite this bracket with the projectors — k) and 
7}(p + k) on the lower and upper mass shells, 



pBC n . 



C x Bp - kBQ 



n-l 



C x Bk 



= ^(p + QBCn-t-'-dBip-k) + ^(p-k)BC n „ 1 ---C x B{p + k) . 
For symmetry reasons, it suffices to show that the first summand of this decomposition, 

{(j> + k)BC n - 1 ---C 1 B(p-k))(x,y) , (3.30) 

is smooth and bounded. 

We proceed in momentum space. We say that a function f(q) has rapid decay for posi- 
tive frequency if it is C , bounded together with its first derivatives (i.e. sup |/|, sup [<9//| < 
oo), and satisfies for every a > the bounds 



sup \uj a f(uj, k)\ 
j>o, feeiR 3 



sup \uj a dif(cu, k)\ < oo . 
j>o, feeiR 3 



(3.31) 



After setting Co = p — k and C n = p + k, the operator product ( 3.30D is of the form 
(1.8). We choose a function g with rapid decay for positive frequency and decompose 
the operator product in the form ( 1 . 12| ) , ( 1.13). It follows by induction that the functions 
Fj all have rapid decay for positive frequency: The induction hypothesis is obvious by 
setting Fq = g. The induction step is to show that for a function -Fj-i with rapid decay 
for positive frequency, the convolution 



dt,)' r dk' 

^- -^Biu-J.k-k^C^iJ^F^iu^k') 



2vr J (2tt) 3 



(3.32) 



also has rapid decay for positive frequency. In Lemma [O], it was shown that Fj is C 1 
and bounded together with its first derivatives. As a consequence, we must only establish 
the bounds ( 3.31 ) for oj > 1. Because of monotonicity u> a < lo^ 3 for a < [3 (and uj > 1), 
it furthermore suffices to show that there are arbitrarily large numbers a satisfying the 
bounds ( p. 31 ); we only consider a = 2n with n G IN. For oj > 1 and u>' G IR, we have the 
inequality 

u 2n < (2cu') 2n e(w') + (2(w - u')) 2n , 

as is immediately verified by checking the three regions to' < 0, < to' < u)/2, and 
u' > uj/2. We combine this inequality with ( 3.32 ) and obtain for u > 1 the estimate 

\u 2n Fj(u,k)\ 



< 



dJ 

W) 

dJ 



dk 1 
(2^)3 

dk' 
(2^ 



B(u -uj',k- k') Cj-xiu', k') \(2uj') 2n @(u') Fj^u', k) 
(2(cu - J)) 2n B{oj -uj',k- k')} Cj-i(u', k') F^u/, k) 



(3.33) 
. (3.34) 



According to the induction hypothesis, the square bracket in ( 3.33j ) is bounded together 
with its first derivatives. Since B has rapid decay at infinity, the square bracket in ( 3.34| ) 
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also has rapid decay at infinity. Thus both integrals in ( |3.33 ) and ( |3.34|) satisfy the 
hypothesis considered in Lemma |Q| for ( |1.10 ), and are therefore bounded. For estimating 
\uj 2n d[Fj\, we differentiate (|3,32|) and obtain similar to (|3.33|) , ( 3.34 ) the inequality 



u 2n diFj{u,k)\ 



< 



+ 



dJ 

dJ_ 

du 



dk' 
(2^)3 

dk' 



(2vr) 3 



diB(u -Lo',k- k') Cj-xiy/, k') [(2u/) 2n 0(o/) F^u', k) 
(2(w - uj')) 2n d t B(u3 — u',k — k')} Cj-ifw', k') Fj-^u/, k) 



This concludes the proof of the induction step. 

We have just shown that for a function g with rapid decay for positive frequency, the 
function 



Fn{q) 



(27T)' 



(BC n - 1 B---BC 1 BC )(q,q 1 )g(q 1 ) 



(3.35) 



has rapid decay for positive frequency. We now consider what this means for our operator 
product (|3.30D in position space. For a given four- vector y = (y°,y), we choose 



g{u),k) = T)(u) e 



-i(uy°-ky) 



where rj is a smooth function with t](uj) = 1 for u < and rj(u) = for to > 1 (this 
choice of g clearly has rapid decay for positive frequency). Since the support of the factor 
Co = (p — k) is the lower mass cone {q 2 > 0, q° < 0}, g(w,k) enters into the integral 
( HID only for negative lj. But for uj < 0, the cutoff function rj is identically one. Thus 
the integral ( [3.35|) is simply a Fourier integral; i.e., with a mixed notation in momentum 
and position space, 

F n (q) = (BC n ^ l B---BC l B( P -k))(q,y) . 

Next, we multiply from the left with the operator (j> + k), 

{{ P + k)BC n „ l B---BC 1 B{p-k)){q,y) = (p + k)(q) F n (q) . (3.36) 

Since F n has rapid decay for positive frequency and (p + k) has its support in the upper 
mass cone {q 2 > 0, q° > 0}, their product decays fast at infinity. More precisely, 

q 1 (p + k)(q) F n {q)\ < const (J) (p + k){q) 

for any multi-index /. As a consequence, the Fourier transform of (3.36) is even finite 
after multiplying with an arbitrary number of factors q, i.e. 



d 4 g 
(2vr) 4 



q 1 (p + k){q)F n {q) 



< const (I) < oo 



for all x and /. This shows that our operator product in position space (|3.3C ) is bounded 
and, for fixed y, a smooth function in x (with derivative bounds which are uniform in y). 
Similarly, one obtains that ( 3.30| ) is, for fixed x, a smooth function in y. We conclude that 
the distribution ( |3.30| ) is a smooth and bounded function. ■ 



We point out that P (x, y) is a non-causal object in the sense that it does not only depend 
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on the potential B in the "diamond" (L^ DiJ) U (L% n L v ) (with and according 
to (2.2)), but on the external potential in the entire space-time. This becomes clear from 
the fact that the support of the operators p(z\, z^) in the perturbation expansion for P he 
is the whole space (^1,^2) £ IP 4 x IP 4 - Especially, it is not possible to express P he (x,y) 
similar to the formulas of the light-cone expansion in terms of the potential and its partial 
derivatives along the convex line xy. 

The non-causal high energy contribution is an effect of higher order perturbation the- 
ory; it vanishes to first order in the external potential 0. According to the decomposition 
into terms of the form (3.30), it comes about because states on the upper and lower mass 
shell are mixed through multiplication with the potential. Qualitatively speaking, this 
mixing only becomes an important effect if the energy (i.e. frequency) of the external po- 
tential is high enough to overcome the energy difference between the states on the upper 
and lower mass shell. This gives the justification for the name "high energy" contribution. 



3.3 The Non-Causal Low Energy Contribution 

In this subsection, we will put the residual argument and the formal light-cone expansion 



of Proposition 3.2 on a satisfying mathematical basis. In order to explain what we precisely 



need to do, we first recall how the light-cone expansion of the Green's functions makes 



mathematical sense: Theorem 2.3 gives a representation of every Feynman diagram of the 
perturbation series ( |2.12 ) as an infinite sum of contributions of the form ( 2.32j) . According 
to the bound (|2.36| ), there are, for any given h, only a finite number of possibilities to choose 
I a and p a ; as a consequence, we get, for fixed h, only a finite number of contributions (|2.32 ). 
Thus we can write the light-cone expansion in the symbolic form 



((-sB) k s)(x,y) = E E ••• sM(x,y) 



(3.37) 



-1 finite 



where '• • •' stands for a configuration of the 7-matrices and nested line integrals in ( |2,32| ). 
According to the explicit formula ( |2.5| ), the higher o-derivatives of S a (x,y) contain more 
factors (y — x) 2 and are thus of higher order on the light cone. This makes it possible to 
understand the infinite sum in (3.37) in terms of Def. |1.2| we can give it a mathematical 
meaning via the approximation by the finite partial sums ( |1.15 ). In Subsection 2J2, it is 
shown that understanding the light-cone expansion via these partial sums even makes it 
possible to explicitly carry out the sum over all Feynman diagrams. 

According to Proposition |3.2| , all the results of Section ^ are, on a formal level, also 
valid for the residual Dirac sea. We begin by considering the light-cone expansion of the 
individual Feynman diagrams in more detail. Similar to ( |3.37 ), the /c th order contribution 
AP rcs to the residual Dirac sea has an expansion of the form 



AP res (x,y) 



E E 

h=— 1 finite 



(3.38) 



where T^ 1 ' is the a-derivative ( 3.14] ) of the lower mass shell T a , ( 3.13j ). In position space, 
T a has the explicit form 



T a (x,y) = lim 



1 



8-7T 3 o<e^o £ 2 - ie£, c 



+ 



327T 3 0<e^0 



lim (log« 2 - ief) + in + cj E 



4? 
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a ^ (-1)3 « 2 y 



with £ = (y — x), c = 2C — log 2 with Euler's constant C, and the function 

n 1 

$(0) = , 3>(n) = 53 r f or n > 1 . 

fc=l k 

The logarithm in ( |3,39| ) is understood in the complex plane which is cut along the positive 
real axis (so that limo< e ^o log(x + is) = log \x\ is real for x > 0). Alternatively, one can 
avoid the complex logarithm using the formula 

lim log« 2 -ie£°) + ivr = log|af 2 | + in e{f) 

0<e^0 

(e is the step function e(x) = 1 for x > and e(x) = — 1 otherwise); thus the complex 
logarithm describes both a logarithmic pole on the light cone and a constant contribution in 
the interior of the light cone. The basic difference between the light-cone expansions ( 3.37| ) 



and ( p.38 ) is related to the logarithmic pole log \a\ in ( |3.39j ). Namely, as a consequence of 



this logarithm, the higher a-derivatives of T a are not of higher order on the light cone. To 
the order 0((y — x) 2 ), for example, one has 



,1 \" ... 1 / ol \ . , , , , 



da) Ta{X ' v) = 32^1^' {a\og\a\) + 0{{y-xY) (,n>2). (3.40) 



This means that the infinite sum in ( 3,3§| ) cannot be understood in terms of Def. 1.2; the 



number of summands is already infinite to a given order on the light cone. In our context 
of an expansion around a = 0, the situation is even worse, because the a-derivatives of 
T a are singular for a — > (as one sees e.g. in (^40|) ). Thus not even the individual 
contributions to the light-cone expansion make mathematical sense. These difficulties 
arising from the logarithm in ( 3.3S| ) are called the logarithmic mass problem (see 0] for 



a more detailed discussion in a slightly different setting). Since we know from Lemma 
P~TT] that the Feynman diagrams are all well-defined, the logarithmic mass problem is not 
a problem of the perturbation expansion, but shows that something is wrong with the 



light-cone expansion of Proposition 3.2 



In order to resolve the logarithmic mass problem, we first "regularize" the formal 
light-cone expansion by taking out the problematic log \ a\ term. By resumming the formal 
light-cone expansion, we then show that the difference between the residual Dirac sea 
and the "regularized" Dirac sea is a smooth function in position space. We introduce the 
notation 

T^( X ,y) = T a (x,y) - log \a\ g^-L^ Li_L (3.41) 

rreg ( = (±) l rf a=Q . (3.42) 



Def. 3.5 The causal contribution p causal to the fermionic projector is obtained from 
the residual Dirac sea P res by replacing all factors in the formal light-cone expansion 
non-causal low energy contribution P le to the fermionic projector is 

given by 

P le {x,y) = P res (x,y) - P causal {x,y) . 
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By the replacement T^ — > T reg ( h \ the formal light-cone expansion of Proposition 3.2 



becomes mathematically meaningful in terms of Def. 1.2. Thus we can restate this result 
as a theorem, leaving out the word "formal." 

Theorem 3.6 (light-cone expansion of the causal contribution) The results of Sec- 
tion [| also apply to the causal contribution to the fermionic projector. More precisely, the 
light-cone expansion of Theorem \2. S\ holds for p causal if we replace by T reg ( h \ ($-4%) , 
and multiply the formulas of the light-cone expansion from the left with the chiral asym- 
metry matrix X . According to Theorem 2.1, all tangential derivatives can be integrated by 
parts. With Def. |/J Def. \2J^, and Theorem \jjUQ, the light- cone expansion can be reduced 
to the phase-free contribution. According to Proposition \2.$ , the phase-free contribution 
consists, to every order 0((y — x) 2g ) on the light cone, of only a finite number of terms. 

We come to the analysis of the low energy contribution. In the following lemma, we 
reformulate the light-cone expansion of Theorem 2.3 in a way where the infinite sums are 
handled more explicitly. 

Lemma 3.7 The light-cone expansion of the k th order contribution ((— s B) k s)(x,y) to 
the perturbation series \2.1i ) can be written as a finite sum of expressions of the form 

OO -1 ny 

E — j F XcCiy-x) 1 [a 1 ,b 1 + n 2 + --- + n k \ m] dz x d\\ □ V% {z x ) 

[ V [a 2 , b 2 + n 3 + --- + n k \n 2 ] dz 2 d% n n z l +q '' V^ C2 {z 2 ) 

[a k , b k | n k ) dz k dit U n * +qk V£> e (z k ) 7 J S^+'"^(x,y) . (3.43) 



ni,...,n fc =0 

X 



' z k-l 

In this formula, we use for the chiral and tensor indices the same notation as in Theorem 
2.5; ai, b\, and qi are non-negative integers. The parameters r and b\ satisfy the bounds 

r < k-\I\ (3.44) 
h > r-l + \I\ , 1 < I < k. (3.45) 



Proof: The form of the expression ( |3.43j ) is straightforward if one keeps track of the 
infinite sums in the inductive construction of Theorem |2.3| ; it is also obvious that we only 
get a finite number of such expressions. The only point which needs an explanation is how 
one can arrange that all infinite sums are of the form 

oo -. 

£-[->-M° n ••• • (3-46) 

For this, we must manipulate the sums when then Laplacian D n is carried out after ( 2.41| ). 
Whenever a Laplacian acts on a factor (y — x) 1 , we shift the summation index by one. 
More precisely, in the case of one factor (y — x), we use the transformation 

oo -. oo -. 

£ - [a, b | n] U n x (y - xf f {n) {x) = £ - [a, 6 | n] ((y - xf U n x f [n) - 2n d'u^f^) 

n=0 n " n=0 n ' 

oo oo 

= (y-^Jl^W,b\n}nlf in) (x) - 2j;- T [a + l J 6 + l|n]nJJ5 i / (n+1) (x) (3.47) 

71=0 71=0 
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(where /( n ) denotes a function depending on n); in the general case of several factors 
(y — x), we inductively apply ( 3.47] ), 

It remains to show that the parameters ai, bi, and q\ are non- negative, and that the 
inequalities ( 3.44 ) and ( [3.45 ) hold. For this, it suffices to consider the leading summand 
Hi = ■ ■ ■ = Wfc = of ( p. 43 ). Since this is a (special) contribution of the form (2.32), we 
can apply Theorem |2.3j with = l{ + rij, bi = ri + n,, and pi = qi. It follows that ai, i>j, 
and (ft are non- negative. For the proof of the inequalities ( 3.44| ) and ( 3.45j ), we proceed 
inductively in the order k of the perturbation theory. For k = 0, we have r = — 1 and 
\I\ = 1, so that the inequalities are satisfied. Assume that (3.44), ( p. 45 ) hold for a given 
k. We go through the construction steps of Theorem |2.3| using the index shift ( |3.47 ) and 
verify that ( |3.44 ) and (3.45) are also valid to (k + l) st order: 

For the proof of ( 3.44 ), we note that additional factors {y — x) are generated at most 
once in the construction; namely, if the derivative <jH x acts on S^' in step 5). The parameter 
r is only increased if either a Laplacian acts on the factor (y 



in step 3) (leading to 

the index shift ( |3.47 )) or if the derivative <f) x does not act on S^ h ' in step 5). In both cases, 
one loses at least one factor (y — x). This gives ([3.44 ). 

For the proof of ( |3.45| ), we take the contribution ( 3.43| ) with n\ = ■ ■ ■ = n/% = and 
apply the construction of Theorem When computing the Laplacian in step 3), we 
shift the index according to ( |3.47 ) whenever a Laplacian acts on a factor (y — x). Denoting 
the number of index shifts by s, we get a finite number of terms of the form 

£ - X c C i$ x S^ +n \x, y) [ V [s, r + s\0]dz(y- x) 1 ^ n n V^\z) 

n=0 n - Jx 

x [\a 1 ,b 1 \0}dz 1 d I z ln? 1 V£l i (z 1 )--- [ V [a fc ,6 fc |0]^^n^^ ) Cfc (z fc ) 7 J . 

Each index shift decreases the number of factors {y — z) and increments the order of the 
mass-derivative of the Green's function, thus 

\I\ < \I\-s , f = r + l + s . (3.48) 

It again suffices to consider the leading summand n = 0; this is a contribution of the form 



( |2.32|) , After extracting the factors (y — x) in step 4), the parameter 6o = r o + n o satisfies 



■Jo 



r + \I + s 



f + \I - 1 



The parameters ty, 1 < I < k, remain unchanged in the construction; they are still the 
same as in ( 3.45| ), 



k > r-l + \I\ > r - (1 + 1) + |J| 



l = l,...,k. 



When the derivative <jD x is carried out in step 5), either r is decremented and |/| increased 
by one, or \I\ is decreased. In steps 6) and 7), the transformations ( p.56| ) and (2.57) may 
only decrease the sum f + |/|. We conclude that, after performing all the construction 
steps, 

k > r-(l + l) + \I\ , l = 0,...,k. 



The index shift /—>/ + ! finally gives the inequalities (3.45) for I = 1, . . . , k + 1. 
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Theorem 3.8 The non-causal low energy contribution P le {x,y) is, to every order in per- 
turbation theory, a smooth function in x and y. 



Proof: We first outline our strategy: According to Def. pTl], P roposition 3.2, and Def. 
the fc th order contribution to P le (x, y) is obtained from fl3.37|) by the replacement S^ h 



AP le (x,y) = X £ £ ... (T^-T 



reg (h) 



(3.49) 



-1 finite 



Because of the logarithmic mass problem, the infinite sum in (3.49) is ill-defined. In order 
to give ( |3.49| ) a mathematical meaning, we manipulate the infinite sum until recovering it 
as a formal Taylor series, which can be carried out explicitly. Finally, we show that the 
mathematical object AP le (x,y) obtained in this way is a smooth function in x and y. 

Consider the light-cone expansion of the Green's functions of Lemma 3.7. Since there 
are only a finite number of contributions of the form ( p.43| ), we can restrict ourselves to one 
of them. In order to get the corresponding contribution to P lc , we replace the factor 
in ( 3.43| ) according to (3.49) by the operator := T^' — T reg ( h \ i.e. more explicitly 



with 



L a (x,y) 



32vH 



lo g|«l E 



{-\y {a{y-xfy 



(3.50) 



We can leave out the factor (y — x) 1 in (3.43) and disregard the chiral asymmetry matrix 
X in (|3.49| ), because they are irrelevant as smooth functions. Furthermore, we can carry 
out the partial derivatives dj, 1 in ( 3.43j ) with the Leibniz rule. According to the chain rule 
dg9p , this may increase the parameters bi\ nevertheless, the inequalities ( |3.44[ ) and ( |3.45| ) 
remain valid. We conclude that it suffices to consider the formal series 



AP le (x,y) 



E 



ni, 



.,n k =0 
ry 



ni! ■ • ■ n k \ 



[ai,b 1 +n 2 + ■■■ + n k \ m] dz 1 D" 1 WAzi) 



■■■ T [a k ,b k \n k ]dz k n n z k k W k (z k )L^ + - +n ^(x,y) (3.51) 
Jzk-i 

together with the bounds ( |3.44 ) and ( ^.45 ), where W\ stand for partial derivatives of the 
potential Vj \ ■ Our task is to give this series a mathematical meaning and to show that 
it is a smooth function in x and y. 

We first consider the case that the potentials Wi are plane waves, 

Wi(x) = e~ iqiX . 

This allows us to explicitly carry out the Laplacians in ( 3.51| ), and we obtain 



AP le (x,y) 



E 



[ai,6i|0] dzi e~ iqiZl 
1 



[o fc ,6 fe |0] dz k e- iqkZk 



Zfc-l 



ni,...,rifc=0 



n k \ 



Ua\ - a-i) pfj 1 ((1 - ai)(a| - a 2 ) 
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l-a 1 )(l-a 2 )---(l-a k . 1 )(a 2 k -a k )p 2 k y k L^+"^(x,y) , (3.52) 
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where ai denote the integration variables of the line integrals (all running from zero to 
one), and where pi are the momenta 

Pi = Qi + (1 - ai) q l+1 + ■ • • + (1 - a t )(l - a i+ i) ■ ■ ■ (1 - a fc _i) q k . (3.53) 

For fixed values of the parameters ai, . . . ,a k , ( |3.52j ) is a product of k formal Taylor series. 
For example by using the continuation of the logarithm into the complex plane 

log \a\ = — lim (log(a + is) + log(a — ie) + 2in) 
2 0<e^0 



(which is, as in ( 3.3S| ), cut along the positive real axis), we can carry out these formal 



Taylor series and obtain 

AP lc (x,y) = [ V [axM^dzxe-^ ■■■ T [a k , b k | 0] dz k e~^ f^-V L a (x, y) (3.54) 
Jx Jz k _ x \daj 

with 

a = (of - ax)p\ + (1 - at)(al - a 2 ) pj 

+ • • • + (1 - ai)(l - a 2 ) • • • (1 - Qfc_i)(a| - a k ) p 2 k (3.55) 

and the momenta pi according to (|3.53j ). This construction is helpful in two ways: all 
infinite sums have disappeared, and the mass parameter a is now in general non-zero, so 
that the a-derivatives of L a are no longer singular. 

After this preparation, we consider the case that the potentials in ( 3.43j ) are 



smooth and satisfy the conditions of Lemma [O] on the decay at infinity. Then the Fourier 
transform V"W is C 2 and has rapid decay at infinity. Since the potentials are partial 
derivatives of V^^ , their Fourier transforms Wi are also C 2 and have rapid decay at infinity. 
The low energy contribution is obtained from (3.54) by integrating over the momenta qi. 
More precisely, 



dbj J (2vr) 4 J (2vr) 4 
[ V [ai, h | 0] dzi e~ iqiZl ■■■ f V [a k , b k \ 0] dz k L a+b (x, y) , b=0 , (3.56) 



where the parameter a depends on the momenta qi via ( |3.55j ) and ( 3.53j ). We must show 



that ( |3.56D is well-defined and depends smoothly on x and y. Qualitatively speaking, 
we can view the (^-integrations as a multiple convolution in the parameter a. These 
convolutions mollify L a +6 i n such a way that the ^-derivatives can be carried out giving a 
smooth function in x and y. Unfortunately, this "mollifying argument" is quite delicate. 
Complications arise from the fact that the (^-integrals are multi-dimensional and that 
we must handle the additional line integrals over a\. The main problem is that the 
dependence of a on the momenta qi becomes singular when the parameters ai approach 
zero or one (see (3.55)). Because of these difficulties, we give the mollifying argument 



in detail. Equivalently to analyzing the regularity in the parameter a, one can take its 
Fourier transform, 

L a (x,y) = / L r (x,y)e" 1 ™ , (3.57) 

and study the decay properties in r for r — > ±oo. We prefer working with the parameter 
r, because this is a bit easier and makes our argument clearer. 



48 



In the first step of the mollifying argument, we transform the (^-integrals into integrals 
over the momenta pi. Since the transformation ( |3.53 ) is volume preserving, we get 



d A Pk 
(2tt)* 



/ [oi,6i | 0] dzx ■■■ ( [a k ,b k | 0] dz k f {a} (pi, . . . ,p k ) L a+b (x,y) | 6=0 (3.58) 

Jx Jzk-i 



with 



f {a} (pi,..., Pk ) := W l {q 1 )---W k {q k )e-^ + - + ^ , 
where p\ and q\ are related to each other via (|3.53j ). The explicit inverse of ( 3,53 ), 

qi = pi - (1 - ai) pi +1 , 1 < I < k and q k = p k , 

shows that this transformation is regular for all values of the parameters a\. As a conse- 
quence, the function f^{pi, ■ ■ ■ ,p k ) is C 2 and has rapid decay at infinity, both uniformly 
in ol\ (more precisely, the Schwartz norms ||/^ a ^||p,2! P ^ 0, are all uniformly bounded in 
ai). Because of this uniformity, we need not care about the dependence of on the 
parameters ol\ in the following. 

According to the bound ( 3.44| ), r < k. Thus we can insert into fl3~58[) the identity 



da r 5{a\ — p\) ■ ■ ■ 5(a r — p 2 .) 



pull the aj-integrals outside and carry out the integrations over p\, . . . ,p r . This gives 



^ ,e <^> - (s)'f 

rv ry 
x [ai,bi | 0] dzx ■■■ / [a k ,b k \0]dz k U(ax, . . . ,a r ,p r+ i, . . . ,p k ) L a+b (x,y) | 6=0 (3.59) 

Jx Jz k _ 1 



d 4 p r 



d A p k 



with 



U 



d A p\ 
(2vr) 4 



d 4 p r 



5{a x -p\) ■■■ 5(a r -p 2 .) f^( Pl ,..., Pk ) 



(3.60) 



In (|3.59[) , the parameter a depends on a/, pi, and a\ via 



ai) a\ + • • • + (1 — ol\) •••(! — a r -i)(a 2 — a r ) a r 



+(l - ai) ■ ■ ■ (1 - a r )(a 2 +1 - a r+1 ) p 2 +1 

+ • ■ ■ + (1 - ai)(l - a 2 ) • • • (1 - a k -i)(a 2 k - a k ) p\ . 

Since is C 2 and has rapid decay at infinity, it follows by evaluating the integrals over 
the mass shells in ( 3.60| ) that U is C 1 in the variables ai and has rapid decay at infinity. 

Next, we take the Fourier transform in the variables a and a\,... ,a r by substituting 
( |3~57l) and 



U(ai, . . . ,a r ,p r+ i, . . . ,pi) 



2tt 



dr r 



2tt 



i >r i Pr+1 1 



,Pk) e 



— i(tioiH \-T r a r ) 
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into ( |3.59| ). The (^-dependence of the resulting expression for AP le has the form of 
plane waves; thus the aj-integrals give ^-distributions. We can then also carry out the 
reintegrations. Finally, the 6-derivatives in ( |3.59 ) give a factor (— zr) r , and we obtain 



AP lc {x,y) = [ V [a 1: bi \0]dzx ■■■ f [o fc , b k | 0] dz k f 

Jx J z k-l J 

r°° dr 



drpr+i f d 4 p k 



(2vr) 4 J (2tt) 4 
x / 7T U(n,. . . ,T r ,p r+ i,. . . ,p k ) (-iT) r L T (x,y) , (3.61) 

where the parameters i\ are given in terms of r and a/ by 

n = (1 - «i) • • • (1 - aj-i) (ai - of) t . 

The Fourier transform L T (x,y) is an integrable function in r which depends smoothly 
on x and y (this can e.g. be verified by writing the infinite sum ( |3.50| ) with the Bessel 
functions J\ and K\ which decay at infinity). Since U is C 2 in the parameters a/ and has 
rapid decay, its Fourier transform U is a function in t\ which decays at infinity at least 
like ©(r" 1 ). 

We now estimate the cty-integrals in fl3.61| ) for I = 1, . . . ,r, from the right to the left. 
Since U decays in r r at infinity like C(r r _1 ), we have the bound 

l 

da r •••U\ Tr =(l-a 1 )..-(l-a r -i)(ar-c4)T ^ SU P ' ' ' (C 1 ~ «l) ' ' ' i 1 ~ Oir-l^y 1 g , 
are[0,l] 

where we have for clarity only written out the a r -integral; g is a function depending on the 
variables n, . . . , r r _i and p r +i, ■ ■ ■ ,Pk- The inequality ( |3.45| ) implies that all factors (1 — 
ai)^ 1 cancel against corresponding factors (1 — a{) in the nested line integrals. The decay 
properties in the remaining parameters T\, .. . ,r r _i remain unchanged by our estimate of 
the a r -integral. Therefore we can proceed inductively in the same way for the integrals 



over a r _i, a r _2, • • • , «i- The bounds ( p.45|) ensure that all factors (1 — ai) 1 drop out. 
Since we get a factors r" 1 in each step, the factor (— zr) r in (3.61) disappears. Using that 



L T is integrable and that we have fast decay at infinity in the variables p r +i, . . . ,p k , all 
the remaining integrals are finite. We conclude that AP le (x,y) is well-defined. 

If we take partial derivatives of (|3,56j ) with respect to x and y, the derivatives act 
either on the exponentials, yielding additional factors ai, (1 — ai), and qi, or they act on 
L a (x,y). Since L a and its Fourier transform L T depend smoothly on x and y, we can 
repeat the above mollifying argument and conclude that AP le (x,y) is smooth. ■ 

For a very rigid mathematician, it might not seem quite satisfying that the light-cone 
expansion of the residual fermionic projector was first performed on a formal level and 
later made rigorous by resumming the formal expansion. We remark that one could 
avoid all formal series in intermediate steps of the construction by already "regularizing" 
the logarithmic mass terms in the Green's functions s after ( |3.12j ). However, this has 
the disadvantage of becoming quite technical. Our procedure is easier to understand 



because we could introduce the residual argument in Subsection 3.1 without entering the 
mathematical details right away. 

Similar to the high energy contribution, the low energy contribution P le (x,y) is non- 
causal in the sense that it depends on the external potential in the entire Minkowski space. 
This can be understood from the fact that the operator products in the perturbation 
expansion ( |3.16|) contain factors p(z\,Z2), whose support is the whole space (z\,Z2) £ 



50 



IR 4 x IR 4 . We point out that, although the statements of Theorem 3A and Theorem ^8 
are very similar, their proofs are completely different. This shows and illustrates that the 
high and low energy contributions describe two different physical effects. In contrast to the 
high energy contribution, the potential in the low energy contribution need not overcome 
an "energy gap;" as a consequence, the low energy contribution plays an important role 
even if the energy of the external potential is small. 

This concludes our analysis of the Dirac sea. We briefly summarize our main results: 



According to Def. |3,1| , Def. |3.3| , and Def. |3.5| , we decompose the fermionic projector in the 
form 

P(x,y) = P causal (*,y) + P le (x,y) + P he (x,y) . (3.62) 
The causal contribution P causal (x,y) has singularities on the light cone, which are com- 



pletely described by the light-cone expansion of Theorem 3.6. The non-causal low and 
high energy contributions P le (x,y) and P he (x,y), on the other hand, are, to every order 



in perturbation theory, smooth functions in x and y (see Theorem 3.4 and Theorem |3,8j ). 

We finally point out that, in contrast to p causal 5 both non-causal contributions P le 
and P he were only studied to every order in perturbation theory, but we did not consider 
the convergence of their perturbation series. For P le , this convergence problem could be 
studied by resumming all phase-free contributions to the formal light-cone expansion of 



Proposition 3.2. However, there seems to be no easy method at the moment to control 
the convergence of the perturbation expansion of the high energy contribution P he (x,y). 



Nevertheless, Theorem |3.4| and Theorem 3.8 give a strong indication that the non-causal 



contribution is a smooth function. Even if singularities or divergences occurred when 
carrying out the sum over the perturbation series (which seems unlikely) , these singularities 
would be of very different nature than the singularities of P causal (x, y). More precisely, 
they could not be expressed in terms of the external potential and its derivatives along 
the convex line xy, because singularities of this type necessarily show up in finite order 
perturbation theory. Thus we can say that, although the non-causal contributions P le and 
phe re q U i re further study in order to get a complete understanding, the decomposition 
( |3.62| ) is well-established. The light-cone expansion gives a method to explicitly calculate 
the causal contribution P causal (x, y). 

A Some Formulas of the Light- Cone Expansion 

In this appendix, we give a compilation of explicit formulas of the light-cone expansion. 
More precisely, we list the phase-free contribution to the light-cone expansion of the 
Green's functions (cf. Def. p.7| ). According to the reduction to the phase- free contribution, 
the light-cone expansion of the Green's functions is immediately obtained by inserting or- 
dered exponentials into the line integrals, see Def. |2.9| and Theorem 2. 10| . Furthermore, 



using Theorem [3^, the formulas can be directly applied to the fermionic projector; they 
then describe the singularities of P(x, y) on the light cone. 

All the following formulas were generated by the C++ program "class_commute" (see 



Subsection 2J3). Our listings are not intended to be in any sense complete; we made a 
selection in order to give the reader a first impression of the form of the singularities. 
The more detailed formulas to higher order on the light cone can be easily obtained with 
"class_commute." Without loss of generality, we restrict ourselves to the left handed 
component of the Green's functions; for the right handed component, the formulas are 
analogous. 
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We begin with the perturbation by a chiral perturbation to first order. The phase-free 
contribution to the order 0((y — x) 2 ) on the light cone is 



phase-free 

XL (s [Xl4r + Xr4l) s )(x,y) ~ 0{{y-x) 



+XL S^(x,y)e / dz [0,1 \0]^A U ) 

J X 

- XL S^(x,y) f dz [0,0 \0}4 L 

J X 

+XL SW(x,y) 4 L (x) 

+ \xl S^\x,y)ij V dz [0,0 | 0] {Ul) 



-XLS<®(x,y)t / dz [1,0 \0](d i A Li ) 

Jx 

+\xl 5 (0) (x, y) i e f x dz [0, | 1] (DA Ll ) 
+XL S^{x,y)e ( V dz [0,1 \l]{^A Li ) 

Jx 

+ XL S^\x,y) dz [0,2 \0](D4 L ) 

Jx 

-2 XL (x, y) [ V dz [0, | 1] {0A Li ) 

J X 

+i xl (x, y) $ jT dz [0, | 1] ($n4 L ) 



- XL S^ l \x,y)i / dz [1,0 \l]{d i UA Li ) 

Jx 

+\ XL (x, y) i e f dz [0, | 2] {U 2 A Ll 



(A.l) 

(A.2) 
(A.3) 
(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A.10) 

(A.ll) 

(A.12) 



where we used the abbreviation £ = (y — x). This formula has the disadvantage that it 
contains ordinary partial derivatives of the chiral potential; it would be better for physical 
applications to work instead with the Yang- Mills field tensor and the Yang- Mills current. 
Therefore, we introduce left and right handed gauge-covariant derivatives D L / R , 



Da 



d 



%A 



Lj 



D 



R 



d 

dxi 



iA 



and define the corresponding field tensor and current as usual by the commutators 



F jk 



3l 



D ck , Ff k 



(c = L or R). 



(A.13) 



In the Abelian case of a single Dirac sea (i.e. / = 1), ([A.13 ) reduces to the familiar 
formulas for the electromagnetic field tensor and current, 



? jk - djA ck -d k A cj , jf = di k A k c -uA c i 



OkAr 



Notice, however, that in the general case of a system of Dirac seas, ( A.13| ) involves 
quadratic and cubic terms in the potential. 

By substituting ( |A.13 ) into ( |A.1| )- (|A.12 ) and manipulating the line integrals with 
partial integrations, one can rewrite the phase- free contribution in a way where the linear 
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terms in the potential are gauge invariant. For example, we can combine ( |A.1| ), ( |A.2j ), and 
( |A.3| ) by transforming the line integrals as 

e F dz [0, 1 | 0] ($A Lk ) = ^ F dz [0, 1 | 0] (y"i$ + d k 4 L ) + O(Al) 

Jx Jx 

= e F dz [0, 1 I 0] 7 J if fc - 4l(x) + I* dz [0, | 0] 4 L + 0(A 2 L ) . (A.14) 

Jx J X 

However, this procedure yields (in the non-Abelian case) quadratic and cubic terms in the 
potential which are not gauge invariant. Fortunately, these gauge-dependent terms are 
all compensated by corresponding contributions of the higher order Feynman diagrams. 
More generally, it turns out that, after summing up the perturbation series for the chiral 
perturbation, we can arrange a gauge invariant phase-free contribution for which the 
insertion rules of Def. 2J} and the statement of Theorem 2.1C hold. This is not astonishing 
in view of the behavior (|2.102 ) of the Green's functions under local gauge transformations; 
we verified it explicitly term by term for all following formulas. We now list the gauge 
invariant phase-free contribution to the order 0((y — x) 2 ) on the light cone. For simplicity, 
we omit all third order terms in the potential which are of the order 0((y — x)°) on the 
light cone and which have a prefactor $ (the combinatorics of the tensor contractions leads 
to many such terms, but they are not very instructive here). 

CO 

xl £((-« (xl4 r + XR4 L )) k v) 

k=0 

phase-free ^ _ ^ + ( ^ _ ^ 

+XL S^(x,y)e F dz [0, 1 | 0] y l F% 

Jx 

+- A XL (x,y)$ £ dz [0, | 0] 7 V Fy 
— XL SW&y)^ 1 ("dz [0,0 | 1] j 

4 Jx 

+XL s^(x, y )e F 'dz [o, 1 1 1] (^; L ; 

Jx 



+XLSM(x,y) I' dz [0,2 |0]j£ 7 fc 



l XLS^{x,y)i f <te[0,0|l](^) 7 * 



2 



— XL S {1) (x, y) $ C J" dz [0, | 2] (pj?) 
-i XL S (0) {x,y)i F dz x [0, 1 | 1] Ffy F dz 2 [0, 1 | 0] F? 

Jx J Zl 

+i XL S^(x, y) Ce F dz, [0, 3 | 0] 7 fc F% F dz 2 [0, | 1] 

J X J Zl 

+i X L (x, y) ee F dz x [0, 2 | 1] j} F dz 2 [0,1\Q} 7 1 

Jx J Zl 

-2i XL S^(x, y) F d Zl [0, 2 | 1] F^ F dz 2 [0, 2 | 0] {$Ff) 

Jx J Zl 

-2i XL S^(x, y) tit? F d Zl [0, 2 | 1] F dz 2 [0, 1 | 0] Ff 

JX J Zl 
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+i XL (x, y) ee F d Zl [0, 2 I 1] 7 fc i^ r dz 2 [0, 2 | 0] 
-\xl S W (x,y)e Fd Zl [0,2|0] 7 J i^ /""das [0,0|0] 7 V^ 
-- XL <S (1) (x, y) C F d Zl [0, 2 | 0] tV i$ r dz 2 [0, 1 | 0] 7 <i^ 

^ J Zl 

+2i X L (Z, y) & F dz, [0, 3 | 0] ^Ff k F dz 2 [0, 1 | 0] F% 

Jx J Zl 

-2i XL (x, y) e F d Zl [0, 1 | 1] f£ F dz 2 [0, 1 | 0] lk F% 

J X J Zl 

~l XL (x, y) $ F dz! [0, 2 I 1] jf F dz 2 [0, | 1] 

^ JX J Zl 

-- XL <S (1) (x, y) $ F dz! [0, 1 | 2] jf F dz 2 [0, 2 | 0] 

^ JX J Zl 

+i XL (x, y) % ee F d Zl [0, 2 I 1] F£ f dz 2 [0, 1 | 1] ) 
+m S« (x, y) jf Ce F d Zl [0, 1 | 2] Ffy F dz 2 [0, 3 I 0] (d k j^) 

F d Zl [o,i 1 2] (a fc ^) f <fe 2 [o, 2 1 o] (d k F l i) 

J X J Zl 



i 

+ 4 
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ry 

/ dzi[0,2 


0]y'7 fc i?* /"^2 [0, | 1] jf 

J Zl 


i 

+ 4 


XL 






ry 

/ dz![0,l 
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l]7 J 7 fe if fc [0, 2 | 0] jf 
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ry 

/ dzi[0,l 


F dz 2 [o, o | o] y 7 fc i$ 

J 21 
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s« 




ry 

/ cfei [0, 1 

>x 


1] 7 V(^) [0,l|0]itf 

J 21 


~2 
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[x,y)U l 


ry 

/ ^i[0,l 


r^ 2 [o,i|o] 7 v(^) 

J 21 



+i XL (x, y) $ Ce F d Zl [0, 1 I 2] (d k rf) F dz 2 [0, 1 I 0] 

JX J Zl 

-iXL (x,y)$ e F dz, [0, 1 | 1] F% F dz 2 [0, | 1] f L 

Jx J Zl 

-iXL S (1) (x,y)$ e F dz, [0, | 2] F% F dz 2 [0, 2 | 0] f L 

Jx J Zl 

+z XL {x,y)$ C F dz, [0, 1 | 1] j{ F dz 2 [0, 1 | 0] F^ 

Jx J Zl 

+2i XL ( x , y) $ ^ F dzi [0, 2 | 1] F& F dz 2 [0, 2 | 0] {d k F l D 

Jx J Zl 

-2i XL {x,y)$ e F d Zl [0, | 2] (d k Ft 3 ) F dz 2 [0, 1 | 0] F kl 

Jx J Zl 

-\xl (x, y) i J" d Zl [0, 1 | 0] y 7 fc if fc £ dz 2 [0, | 0] 7 ' 7 m F^ 
+3i X L y) i F dz 1 [0, 1 | 1] F& F dz 2 [0, 1 | 0] F kl 

Jx J Zl 
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Ik 



dz 2 [0, 1 | 1] F, 



L 

mj 



dz 3 [0, 1 | 0] FT 



-2XLSW(x,y)ti?Z k / ^![0,4|0]yF ; 

Jx 

-2 XL S^(x, y) U j i k F dz! [0, 3 | 1] F dz 2 [0, 3 | 0] F^ F dz 3 [0, 1 | 0] F, 

Jx J Z\ J Z2 

-2 X l S^(x, y) U j t [ V dzx [0, 3 | 1] F^ k F dz 2 [0, 3 | 0] 7 l F^ F dz 3 [0, 1)0]^ 

Jx J 21 J Zo 



rnt 
L 



nit 
L 



-2 XL (x, y) F dz x [0, 3 | 1] F^ k F dz 2 [0, 3 | 0] F? j F dz 3 [0, 1 | 0] ^itf 

Jx J Zl J Z2 

We call this formulation of the phase-free contributions purely in terms of the Yang- 
Mills field tensor and the Yang- Mills current the gauge invariant form of the light-cone 
expansion. 

For clarity, we mention a subtlety in the transformation to the gauge invariant form: 
The gauge invariant phase- free contribution implicitly contains tangential derivatives of the 
chiral potential (as one sees by writing out F c and j c in in terms of the chiral potentials). 
Thus it is not a phase- free contribution in correspondence with Def. 2/7; as a consequence, 
it is not obvious that the statement of Theorem 2. 10| holds. In other words, one must 
be very careful when transforming the line integrals in order to ensure that the insertion 
rules of Def. 2.9 and Theorem 2.10 remain valid. A safe method is to insert all ordered 
exponentials into the line integrals before performing the partial integrations. According 
to Theorem p. 10 , the phase-inserted formulas coincide with the light-cone expansion of 
the Green's functions. Therefore the partial integrations become identical transformations 
of the Green's functions; we need not worry about the insertion rules. In the final step, 
we again take out the ordered exponentials from the line integrals, verifying that they are 
still in accordance with Def. In the example of the line integral ( A. 14 ), the partial 
integration of the phase-inserted contribution gives 

e fc F dz [0, 1 I 0] Pe -< X>£ @A Lk (z) Pe^SH^ ^~ z ^ 

Jx 

= e F dz [0, 1 | 0] Pe-C a l (>?Ff h + i [4 L , A Lk ] + d k 4 L ) , Pe"*/" A ™ {y ~ zU 

Jx ^ ' \ z 

dz [0, 1 | 0] Pe^r^ (*"*>« (>F/ fc + i [4 L , A Lk \) z P e ~ l H A L (v-*U (a.15) 



*T (y-z)r, 



-4 L (x) Pe-*/"^-*)" 
+ F dz [0,0 | 0] Pe-*£^ {z - x)l 4 L {z) Pe~ l l 

Jx 

H k F dz [°> 1 I °] Pe_l/: A ' L {Z ~ Xh i i A Lk, 4l) (*) Pe" 4 ^ 

Jx 



(A.16) 
(A.17) 

(A.18) 



Thus the correct transformation of the phase-free contribution is 



tt / dz [0, 1 | 0] $A Lk ) 



phase- free 



t I dz [0, 1 | 0] 7*i$ + / dz [0, | 0] 4 L - 4 L (x 



Notice that ( [A. 18 ) are the derivative terms of the ordered exponentials; they get lost if 
the phase-free contribution is transformed in a naive way. 

It remains to consider the scalar /pseudoscalar perturbation; i.e., we must study how 
the dynamic mass matrices Y^m(x) show up in the light-cone expansion. We begin with 
the case of a single mass matrix. To first order in the external potential, the corresponding 



55 



Feynman diagram has the light-cone expansion 



phase- free mf/ \ 2 \ 



XLm(-s(-XLY R - xrXl) s)(x,y) x 0((y - x 
+\xLm sW(x, y) $ £ dz [0, | 0] ($Y L ) 

+XL m S<®(x,y)Y L (x) 
+ XL mS^(x,y) [ V dz [0,1 | 0] (DY L ) 

J X 

+\ XL m S (1) (x, y) jf £ dz [0,0 \1]($C\Y L ) . (A.19) 

Similar to ( |A.l| )-( [A.12j ), this formula involves partial derivatives of the potential, which is 
not a gauge invariant formulation. Since the chirality changes at every mass matrix (see 
e.g. Def. |2.9[ ), the correct way to make the light-cone expansion gauge invariant is to work 
with the gauge-covariant mass derivatives D, A given by 

(DiY L ) = Df Y L - Y L Df = (diY L ) - iA Li Y L + %Y L A Ri 
(AY L ) = D^D'Yl) = D\ (D'Yl) - (D'Yl) Df , 

and similarly for Y R and higher derivatives. Rewriting ( A.1S| ) with gauge-covariant mass 
derivatives yields additional terms which are linear or quadratic in the chiral potentials 
and which are not gauge invariant. But, similar as described for the chiral perturbation 
above, all these terms cancel if the sum over the perturbation series for the chiral potentials 
is carried out. To the order 0((y — x) 2 ) on the light cone, we obtain in this way the gauge 
invariant phase-free contribution 

oo 

XL™ ((- s (xL4 R + XR4 L )) ni s(xLY R + XRY L )s((-XL4 R -XR4L) s ) n2 )( x ^y) 

Hi ,712=0 

= XL £ ((-s4 L rsY L s(-4 R sr)(x,y) Pha x free 0{{y-xf) 

711,712=0 

+ \ XLm S-(°) (x, y) % [ V dz [0, | 0] V (DjY L ) 
+XL m S(°Xx,y)Y L (x) 
+XLmSV>(x,y) f dz [0,1 \0] (AY L ) 

J X 

+\ XL m (x, y) i £ dz [0, | 1] 7 J (D 3 AY L ) 
+ % -XLm (x, y) $ C [ V dzi [0, 1 | 1] j[ f dz 2 [0, | 0] 7* {DjY L ) 



-i XL m SW(x, y) i C f dzi [0, 1 | 1] F$ f dz 2 [0, 1 | 0] 7fc (D k ^Y L ) 

J X J Zl 

"7 XL m (x, y) i f V dz, [0, 1 | 0] 7*7^ f dz 2 [0, | 0] 7 J (DjY L ) 



-i XL m SU(X, y) C I" dz x [0, 2 | 0] 7 fc i^ f dz 2 [0, | 0] y (D^) 

J X J Zl 

-\xLm ( X) y ) £ dz [o, 1 | 0] 7 V'i$ Y L 
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"2 XL m 



sW(x,y)$ F dz [0,0 \l]j^ 7 k Y L 

J X 

-i X L m S<U(x, y) $ J V dz [0, | 1] ^F§(D*Y L ) 
+ l - xl m SV> (x, y) jf ? jf" ttei [0, 2 | 0] f (Z^) jT dz 2 [0, ] 1] j 
+^ XL m (x, y) ({ f f d*i [0, 1 | 1] 7* /" tfoa [0, 2 | 0] if 

* Jx J Z\ 

-i X L m y) £ r f dzi [0, 1 | 1] lk (D k ^Y L ) F dz 2 [0, 1 | 0] F$ 

Jx J Z\ 

-I- xl m y) $ F dz! [0, 1 | 0] 7 J (D,Y L ) F dz 2 [0, | 0] 7 fc 7 ^ 

+i XL m S^(x, y) C F dz, [0, 2 | 0] 7* (D,^) F ttea [0, 1 | 0] 7 fc F* 
-2i XL m 5«(x, y) f T [0, 2 | 0] A" dz 2 [0, 1\0]F, 

Jx J Zl 

+i XL m SF> (x, y) $ F dz [0, | 1] (D j Y L ) fF, 

J X 

% xl m ( X) y ) $ F dz [0, | 1] Y L jh k 

Jx 

sW(x,y)?Y L F dz [0,0 | 1] 

Jx 



■ ji 



2 

-iXL m 



1 xl m (x, y) F dz [0, 1 | 0] Y L y 7 *F$ 



2 



- ~ XL m ( x , y ) Y L jT [0, | 0] 7Vi$ 

+XL m S^(x, y) i F dz, [0, 2 | 1] F dz 2 [0, 2 | 0] Ff f [0, | 0] 7 l (£>,YL) 
+XL m S«Oz, y) £ e J T dz! [0, 1 | 1] F%. F dz 2 [0, 1 | 0] lk F? k Y L 

Jx J Zl 

+XLmS^{x,y)U i e Fdz! [0,2|1]F£ F dz 2 [0, 2 | 0] V (D t Y L ) F[0,l\0]F ki 

Jx J Zl J Z2 



-XL m 



^ S {1) (x,y)f C F dzi [0, 1 | 1] Ft F dz 2 [0, 1 | 0] F L 7 ,- F« 
+XL m SW(x, y) i & f tfax [0, 1 | 1] F L f dz 2 [0, 1 | 0] F$ 

Jx J Zl 

+XLmSW(x,y)Uit; j F dz, [0, 3 | 0] 7 ; (D t Y L ) F dz 2 [0, 1 | 1] f]g /"[0, 1 I 0] i$ 
+XL m ^(a, y) £ Cif /" tfax [0, 2 | 1] 7 ; (A^l) /" ^2 [0, 3 | 0] F* F[0, 1 | 0] F k R 

Jx J Zl J 22 

+XL m S^(x, y) % e< f <tei [0, 1 I 1] *l /" ctea [0, 1 | 0] 

Jx J Zl 

-2 XL m ( x , y ) tfJ Y L F dz, [0, 1 | 1] F* F dz 2 [0, 1 | 0] F% . 

Jx J Zl 



Zl 

2> 



The higher orders in the mass matrices are treated similarly. To the order G((y — x t . 
on the light cone, only the terms up to fourth order in m Y^m contribute (see ( |2.95| )). One 
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gets the formulas 

CO 

XL m 2 £ ((" s 4 L ) ni s Y L s (-4 R sr Y R s (-4 R sD{x,y) 

ni,n2,n3=0 

phase-free * m 2 g (0) (x>y) j [" dz %0 \ 0]Y L Y R 
* Jx 

-^XLm 2 S^(x,y)i Fdz 1 [0,l\0]^(D j Y L ) F dz 2 [0, | 0] -y k (D k Y R ) 

^ Jx J z\ 

+\x L m 2 S« (x, y) $ J" dz [0, | 1] Y L (AY R ) 
+i XL m 2 ( x , y ) % F dz [0, | 1] (DiY L ) (D' l Y R ) 

Jx 

+\x L m 2 SW (x, y) $ J V dz [0, | 1] (AY L ) Y R 
+i X L m 2 ( Xj y ) F dz [0, 1 | 0] Y L ^(DjYr) 

Jx 

+i XL m 2 (x, y) F dz [0, 1 | 0] j j (DjY L ) Y R 

Jx 

-iXL m 2 ( Xj y ) Y L F dz [0, | 0] 7 J (DjY R ) 

Jx 

-\x L m 2 sW(x, y) % e F d Zl [0, 1 | 1] $ F dz 2 [0, | 0] Y L Y R 

Jx J Z\ 

+ XL m 2 sW(x, y) $ C F dz! [0, 1 | 1] F± F dz 2 [0, 1 | 0] Y L (D^Y R ) 

Jx J Z\ 

+XL m 2 S (1 \x, y) $ e F dz, [0, 1 | 1] F± F dz 2 [0, 1 | 0] {LPY L ) Y R 

Jx J Z\ 

+\ XL m 2 (x, y) 4 jT dz, [0, 1 | 0] f dz 2 [0, | 0] Y L Y R 

+XL m 2 S^(x, y) C F dz x [0, 2 | 0] F dz 2 [0, | 0] Y L Y R 

Jx J Z\ 

-\xLm 2 (x, y) 4 C F d Zl [0, 2 | 0] Y L Y R F dz 2 [0, | 1] jf 

£ Jx J Z\ 

-\xLm 2 S& (x, y) 4 C F d Zl [0, 1 | 1] Y L Y R F dz 2 [0, 2 | 0] jf 

£ Jx J Z\ 

+XL m 2 S^(x, y) 4 e F dz! [0, 1 | 1] Y L (D^Y R ) F dz 2 [0, 1 | 0] i$ 

Jx J Zl 

+ XL m 2 SW(X, y) 4 C F dz, [0, 1 | 1] (D^Y L ) Y R F dz 2 [0, 1 | 0] F± 

Jx J Z\ 

+- A xl m 2 sW(x, y) 4 f x dz! [0, 1 | 0] Y L Y R J" dz 2 [0, | 0] 

+XL m 2 S^(x, y) C F dz, [0, 2 | 0] Y L Y R F dz 2 [0, 1 | 0] + 0((y - x) 2 ) 

Jx J Zl 

00 

XL m 3 J2 i(~ s 4 L ) ni * Y L s (-4 r sr Y R s (-4 L a y» Y L s (~4 R sD(x, y) 

ni,ri2,n3,n4=0 
phase-free , <,« fV 



XL m 3 (x, y) Y L F dz [0, I 0] Y R Y L 

Jx 
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+ 1 XL m 3 y) % f V dz x [0, 1 | 0] Y L Y R dz 2 [0, | 0] j j (DjY L ) 

J X J 21 

+^ XL m 3 5«(a;, y) % dz x [0, 1 | 0] ^(DjYl) jf" dz 2 [0, | 0] Yr 17, + 0((</ - x) 2 ) 

oo 

ni,...,n 5 =0 

x Y R s (-4 L s) n * Y L s (-4 R sr Y R s (-4 l s) n5 )(x, y) 

phase-free % , , , , fV , , /"V , , ,, l0l 

- 5 XL m 4 S« ( z , y) ({ / [0, 1\0]Y L Y R dz 2 [0, | 0] Y L Y R + - x) 2 ) 

^ J X J Z\ 

These classes of Feynman diagrams completely characterize the Green's functions to the 
order 0((y— x) 2 ) on the light cone; notice that we only get a finite number of contributions. 
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